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CHAPTER I 
INTRODUCTION 

Knowledge of the space environment is a necessity for the design and 
operation of space behicles, the space shuttle/ and Spacelob. Complex physical 
characteristics of the space environment are essential in providing a systematic 
background for the future advancement of knowledge* 

In this report, we shall present the results of our studies during the 
period March (970 ~ February 1974, of this contract. We have covered a wide 
range of problems in the space environment, such as the problems of the dynamical 
behavior of the thermosphere, hydromognetic wave propagation in the ionosphere, 
and interplanetary space environment. The theories we used to analyze these 
problems range from a continuum theory of magnetohydrodynamics to the kinetic 
theory of free molecular flow. This is because the problems which we encountered 
cover the entire range of the Knudsen number (i .e., the ratio of mean free path 
to the characteristic length). From these studies, significant results were obtained. 
They are summarized as follows. 

In the study of the dynamical behavior of the thermosphere, a theoretical 
model based on magnetohydrodynamic theory is suggested- We have shown that 
it is possible to predict the dynamical responses of the thermosphere due to a 
geomagnetic storm by using this model . It was demonstrated in this study that 
the quantities of wind speed, electric field strength, and joule heating due to a 
geomagnetic storm, which we predicted, agreed well with the observational data 
obtained by LOGACS (Low-G Accelerometer Calibration System) flown on a 
Agena satellite of the United States Air Force. Details of this investigation are 
given in Chopter II. 

During the storm time, the ionospheric heating is a very important problem, 
which may trigger hydromognetic waves. However, study of the propagation of 
hydromognetic waves in the region of the ionosphere is a very complex problem due 



to the vanations of the collision frequencies of ions# electrons, and neutrals. 
Therefore, we shall restrict ourselves to dealing only with the propagation of 
hydromognetic waves in the upper F2 region of the ionosphere. In this region, 
we find that electrons are in a transitional regime from collisional to 
coll isionless conditions and ions are in a collisionless regime. Thus, we have 
derived a set of governing equations which are based on the fact that the 
isotropic electrons are fluid like, and the anisotropic ions follow kinetic 
equations modified by lon^electron collisions. It has been shown that the 
magneto-acoustic waves of a period of ^ 10 seconds are dissipated by ion Landau 
damping and electron thermal conduction and viscosity. The numerical solutions 
under ionospheric conditions show that electron heating has always prevailed 
over ion heating in agreement with observations in the present study. The 
implications of possible mirror and/or firehose instabilities, which may cause 
turbulent or Irregular fluctuations of fields, are discussed. The details of this 
study are included in Chapter III . 

In Chapter IV we investigated the properties of the interplanetary medium, 
particularly that of the solar wind interaction with planetary objects. Because of 
the physical nature of the interplanetary medium (i .e., very large Kundsen 
number) a purely kinetic treatment is used for this study concerning the inter- 
action of the solar wind with any "small" planetary object. "Small" refers to the 
coses where the thermal gyroradius of protons In the solar wind Is arbitrarily 
taken to be > O.I times the radius of the object under investigation. The "object" 
may possibly include an ionosphere or magnetosphere- The col lislonless Boltzmann 
equation, neglecting the magnetic field, is used to calculate steody-stote 
profil es of density and velocity around the object. In order to include the effects 
of surface reflection from the object, a finite, angularly distributed potential 
field on the surface of the object Is also considered. A symmetric potential field 
is also examined for the purpose of comparison. A low density plasma void in the 
umbra! region and a compression (above the value of the upstream undisturbed density) 
in the penumbral region are clearly found. The present technique, despite Itsneglect 
of the interplanetary magnetic field. Is proposed as an alternative zeroth order 
approach to the continuum, local magnetic anomaly, and griding center approoches 


2 



used by others for the particular case of the moon. Some recent, potentially relevant 
observations on, and in front of, the moon are discussed. Of particular interest 
is the development of a high density gradient which emanates from the limit and is 
inclined in the direction of what is known in continuum theory as the "Mach cone.” 
This structure leads to the following speculation: if the object’s size were to increase 
(or, alternately, the proton gyroradius were to decrease) by one or two orders of 
magnitude, a shock wave would develop first at the limb, gradually building up 
to a conventional detached bow shock . Thus, gasdynamic theory ~ as we think of 
it in the Eorth's case - can be constructed as required, on the basis of kinetic theory. 

Finally, we shall outline a recommendation for future study in Chapter V - 
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CHAPTER II 

DYNAMICS OF THERMOSPHERE 


II- 1 Introduction 

A pnomagnetic storm is a natural disturbance of the geomagnetic field on ond 
above the earth's surface. Its origin and associated phenomena are very complicated. 
Theoretically, a variation of the magnetic field may be caused by the following phenome- 
na 1) variation of the geomognetic poles, 2) presence of a local inhomogeneous 

^ "b 

electric field -r — ~ ^ x E , and 3) variation of electric currents. 


The magnetic field generated by the earth ‘s magnetic poles is relatively steady. 
Its variation Is negligible in the period of a geomagnetic storm. On the surface of 
the earth/the effect of the Ionization of the atmosphere is small and there is no signifi~ 
cant electrical field that can be detected. Therefore, it is believed that the geomag- 
netic disturbance recorded on the earth's surface is mainly due to the existence and 
varialion of electric current flowing obove the earth (Chapman and Bartel f 1 J ). 
Indeed, if the electric current data were given for the whole environment of the 
earth, then the magnetic field could be calculated everywhere, at least theoretically. 
This would explain the disturbance of the magnetic field recorded on the earth's surface 
during the geomagnetic storm. Actually, there are always electric currents flowing 
above the earth, and their intensities and directions are changing from time to time. As 
a consequence, the geomagnetic field values recorded on the earth will never be 
constant. If the transient geomagnetic field variations are smooth and regular, we say 
it is quiet; otherwise, It is said to be disturbed. 

On a quiet day, the magnetic voriatlon proceeds mainly according to local 
solar time with a small part related to the moon. The two parts are called solar daily 
and lunar daily magnetic variations, and their correspondence to magnetic fields is 
denoted by Sq and L. During a mognetic disturbance, additional electric currents 
flow in the ionosphere. They are superimposed on the Sq and L currents , From the 
analysis of geomagnetic disturbance data, we see that at least five components of 
electric currents are involved in a magnetic storm. They are described as follows: 
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1 • DCF * disturbance due to current attributobie to solor corpuscular flux 
at magnetopause, 

2. DR “disturbance due to magnetic ring current, 

3* DP “polar current, which is strongest in auroral electrojets, 

4. DT “ disturbance due to magnetospheric tail currents, ond 

5. DG “ disturbance due to induced ground currents. 

Since the electric currents ore flowing in the slightly ionized gas medium of 
the upper atmosphere, the other ossociated physical and dynamic effects such os wind 
generation, joule heating, and density and temperature variations can be predicted 
by mognetogasdynamic theory (Piddington [2 ])- A theoretical feature of the geo- 
magnetic storm and its associated phenomena are illustrated in Diagram I. (See also 
[2], P. 13.) 

Therefore, o complete analysis of the problem of the dynamic structure of the 
upper atmosphere due to a magnetic storm has to Include consideration of the distur- 
bance currents and their associated electric and magnetic fields and the dynamic quali- 
ties simultaneously since they are interrelated with each other (see DiagramI). However, 
the solution of this problem is difficult due to the complexity of the mathematics. Many 
authors have studied one or a few parti culor effects separately and made assumptions 
about the other physical qualities rather arbitrarily; Cole SJstudied the joule heating 
of moving ionized gas on the assumption of steady and uniform electric and magnetic 

fields* He concluded that the joule heating effect is significant during a geomagnetic 

-5 -3 ~ 1 

storm. The joule heating may be on the order of 10 erg cm sec in the region of 

5 

100 to 200 km, and the wind as high as 10 cm/sec, if the assumed magnetic and elec- 
tric field configurations are correct. Later, he [4] extended the results by Including 
viscous effects. Maeda and Koto [5 jhove given an excellent review of the problems 
of electrodynamics of the ionosphere in which the problems of conductivity, wind and 
the dynamo theory, drift and its effect on the ionospheric formation, and the interaction 
between wind and electromagnetic field ore discussed in detail . Thomas [6], and 
Thomas and Ching [7 ], opplying a one-dimensional vertical model, reproduced the 
height profile and the mean time lag of the density disturbance by assuming; that the 
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heat input due to a magnetic storm is given* VoHond and Mayer [8 ] reanalyzed the 
same problem using a three-dimensional thermospheric model. 

In this study, a theoretical model is established for calculoting the joule heat 
ing and winds from the geomagnetic variations recorded at storm time. The mathe- 
matical formulation is based on magnetohydrodynomic theory (Chang, Wu, and 
Smith [9] ). Foraday's law is employed for determining the electric field from the 
magnetic field disturbonce data during the geomagnetic storm. 
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11^2 Basics of Theoretical Model 


; 


II-2-1 Basic Equations 

The thermosphere may be considered as o oont^oous medium with a finite 
electric conductivity (Cowling [10]). The presence of electric and magnetic fields 
in the conducting medium will give rise to two principal effects; First, body force 
(lorentz force) and, second, energy generation (Joule heating). These must be taken 
into consideration in the momentum and energy equations. We will derive the equations 
of mass, momentum, and energy conservation as follows: 

Mass Conservation 

The equation of mass conservation is the same os in ordinary fluid dynamics, 

namely 

-4^ + V • ( p V ) = 0 (1) 

at 


where p denotes the mass density and v the velocity vector. 

Momentum Conservation 

The equation of motion of a continuum medium in general can be written 
as (Couchy equation) 


P 



D t 


?iP. . 

——Li — + p 
ax. ^ 
j 


( 2 ) 


where v. are the components of the velocity vector v, P.. are components 
I * J 

of the stress sensor, and F. denotes the components of the body force F. For a 
Newtonian fluid the stress tensor can be expressed os follows; 


P.. 

'J 




( 



ax. 

j 


■t 


Sv. 



(3) 


where p is the static pressure and is the viscosity. By substituting Eq. (2 ) 
into Eq . ( 3 ) we obtain 



V p + X + F 


(4) 
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where 


X = - ’ ( Ti y • V ) + n V + 7 (v • V ) ] 


+ 2[((yn) • y)v + (vTi)x(vxv)] 


and fhe body force F Is 


F -p E+jxB+pg 
6 

— # — + 

where p is the charge density, E is the electric field strength, B is the 
e 

magnetic induction, and j is the sum of the conduction current and the current 
flow due to convective transport of charges » 

Energy Conservation 

The rate of increase of total energy irj the fluid of a moving volume O 
is given by 


f_Di 


d a 


r _L sieSd a +r 

J 2 D t j ' 


d a 


J 

where e is the internal energy per unit mass, 

In order for energy to be conserved, this must equal the energy inputs per 
unit of time from other sources. These are 
1 . Joule heoting ” J (6 • j ) d a 


-f V’ (XVT )da, where X 
_ /* V’ 

done by surface force — - # P. . d s, applying Gouss* 

^ ' U J 


2. Heat conduction per unit time 
is the thermal conductivity. 

3 . W ork 
theorem 


- f Z V. P.. ds. = - r Z ^ (v. P.. )d a 

J i,j 1 IJ J J .. ax , ,J 
S V 


The energy equation is obtained as 


r a = / E -7da-/ 7- ( 

V w \/ 


J 


X 7 T ) d 0 


V 


. . a y ( V. ^ 

‘ ‘J 


(5) 
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The last summation on the right hand side of £q. (5) can be rewritten as 

-Z^(pv,).cp 

I I 

where, from £q. (3) 



Maxwell Equations and their Approximation 

The Maxwell equations are needed to determine the electromagnetic 
quantities* Thus, 


7 • B = 

0 


(6) 

7 • i ^ 

1 

G e 


(7) 

7 X E = 

hi 

bt 


(8) 

7x B = 

M ( J + 8 
0 

at ^ 

(9) 


and Ohm's Law; J“<^(H+vx B)+ 02 (E+ v x B)j/| B | (10) 

A number of approximations can be made which will .be valid as Tar as the present 
analysis is concerned: 
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^ E 

1) The displacement current e may be neglected in Eq«_(9). In order 

S E 

to show that valid approximation/ we shall compare e — to the conduc- 

tion current J = a Eo If E Is assumed to vary periodically with time with frequency UJ , 

BE 

then the ratio of the amplitude of e /B t to the amplitude of J is approximately 

•(-If), 


g(BE/Bt) max _ 
( o E ) max 


e 0) 




fmax 


max 


eU) 


-3 


The value for O in our case is about 10 mhos/m» The value of s in a vacuum 

A 2 “2 -1 

is approximately 9 X 10 forads/m. is about 10 sec , This gives 


e ( B E/^t) 


max 


= 10 


-10 


( a E )max 

which shows that the displocement current can be indeed neglected. 

2) The electrostatic body force can be neglected* 

If E varies Iineorly over some small region, then from Eq. (7) we see 

e E 


Therefor© 


" L 


i- ^ 

p E = 
a 


e £ 


In comparison to the Lorentzion force JxB^o^vB ,we have 

2 


e E 


P E 
e 


^ ® at V 


2 „2 

ev B 

ol V 


e v 

oT 


The values in our problem are estimated as follows 

e - 10 farads/m 
V = 10^ m/sec 

L - 10 m 
~ -3 

cr == 10 mhos/m 


® V ~ 


L a 


10 


-n 


Therefore, the electrostatic body force can be neglected in comparison to the Lorentz 
force. 
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Now, let us recollect the above derived equations os follows: 
Maxwell 's equations: 


V X E 


d B 
^ t 


( 11 ) 


V X B = M j 
o 

7 - 6^0 
e V • E = P 

e 


02 ) 

( 13 ) 

( 14 ) 


Conservation of Moss: 


-|-f- + 7 • ( p ; ) = 0 (15) 

Equatioi^ of Motion: 

P -§^ = -7p+ x+jxB (16) 

Conservation of Energy: 

D $ -* -* 

P -57- = E • j + 7 (X 7 T ) -7- (p^ ) + tp ( 17 ) 


Ohm's Low: 



E’ 


+ a 


•*4 *4 

B X E’ 

"IW 


( 18 ) 


where ^ ' 


E' = E + 


X B. 


This set of equations will form the basis of our theoret- 


ical crnalysiSe Various opproximations will be introduced. 
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11 - 2^2 Som© Fundamental Solutions 

S!ne@ the mggrsetohydrodynomic equaHons [£qs o (11) - (18) ] combine 

the full complexity of Maxwell's equations and the fluid dynamic equations, it is 
obvious that they will be extremely difficult to solve in the general form. Exact 
solutions exist only for a few special cases. We will discuss two of them; namely, 
the Piddington [11] solution for a non-viscous gas and the Hartman solution for a 
viscous fluid [12]» They do not necessarily correspond to the conditions during mag> 

netic storm; however, the solutions will give us some insight into the physical nature 
of the problem. 

Consider 0 uniform ionized gas in steady and homogeneous electric 
and magnetic fields. They are assumed to have the following forms 


E = E ( S , E , 0 ) and B = ( 0, 0, B ) 

X y z 

We further assume that the pressure gradient is 


l£_ 


where I denotes the unit vector in the 


x-di recti on. 

By applying Eq. (18) and Eq. (16) and neglecting the viscous term, we 

obtain 

— * 

B X E' 


j == a. E* a. 


2 B 


(19) 


and 


D ] X B - 7p 


( 20 ) 


Dv 


At hydrostatic equilibrium -gy ® ^ 


T ^ ^ P X B 

^ 

z 


= a, E- . 


B X E' 
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It follows that 


E / B - 


X X z 



- E /b + 3xcr , 

X Z 2 /n ^ 

^ ' B a, a 

z 1 o 

( 21 ) 


J = 0 , 

X 


= ^ p 


h X 


( 22 ) 


whore tj = ^ '^l 


Therefore/ the current j and velocity v are obtained in terms of E/ B and Vp . The 
Joule heating is given by 


Q = J • E' == j • (E+ V X B) 


(23) 


This meons that if we simulate the mechanism of geomagnetic activity by switching 
on the electric field E / after a period of time the velocity v and Joule heating Q 
will reach the volues given by Eq» (21) and Eq« (23)* 

(B) Hartmann Solution 

Cole [4 ] studied the problem of heating and dynamics near auroral electro- 
jets by considering the following model as shown below: 

aZ 

— * 

B 

■ - ” ' * 

-L L X 

The directions of the electric and magnetic fields are shown in the figure and the direction 
of flow is perpendicular to the x - z plane* If the viscous effect of the fluid is taken 
into consideration/ the momentum equation can be written as: 



V 

X 

B x^ 




- j B = 0 
■'x 


(24) 
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where j (E + v B ) » 

■'x 1 ' y 

Equation (24) can be rewritten as 


d V 


X 


.2 
d V 


E||B 


ax 


2 


(25) 


Let us consider the thickness of the sheet of the Qurorol zone to be 2L and 
assume that the motion of the gas outside this sheet is smalL We can then set the 
boundary condition as v sa 0 at x ~ ± L, Then the solution of Bq. (25) that sati 
fies this boundary condition is given by 


where 

as follows: 


V (x) 

y 


111 


1 - 


cosh (M| x/L)’ 


cosh 


(26) 




t (^1^^ ) B The heating within the sheet can be expressed 


Q 


^ 2 cosh (M^ xA ) 
cosh M 


(27) 


1 


n~2~3 Time Dependent Problem 

In the last section we have discussed steady solutions in mognetogosdynamlcs 
which indeed con help to understand the baste mechanism of Joule heating during a 
geomagnetic storm^ However^ they ore not sufficient to explain the time variation of 
the geomagnetic disturbance^ Actually all the quantities -electric field, magnetic 
fiel4ond velocity-are time dependent during the entire period of a geomagnetic storm. 
Therefor©, the analysis for the unsteady case is necessary. 

Let us recoil the equation of motion, namely Eq« (16) 

-4 

n w -4 -4 

. -4 

Suppose that the current data jare known (by either observation or theory). Thofi 
the induced magnetic field may be calculated by £q. (Q), together with the earth's 
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magnef-Jc field. Therefore/ the lost term of (16) becomes known. The problem 
is then reduced to a dynamic problem of the atmosphere with on additional given 
driving force (Lorentz force) and heat source (Joule heating). 

On the other hond, if the magnetic field data in the thermosphere are available, 
we will simplify the problem by eliminating the variables E and j from the system 
of equotions* This is done by solving Eq. (18)/ namely/ 

E = vxB+<x(7xB). (28) 

where a = ^ 2 ^^^ I ^ 

Substitution of j from Eq. (12) then gives 


V X B 


- V X B + 0^ 


Vx B 


(29) 


which/ when combined with Eq* (11)/ gives 

+ 9 X - 9x (v X B) + V X — (Vx B)x B = 0 . (30) 

o t [JO |i 

o I o 

We notice that Eq. (30) Includes only two variables B and v . If B (x, y, z, t) 
is assumed to be given/ then the velocity vector can be solved from Eq. (30) with 
properly imposed boundary conditions. 

Equations (11~18) may be reduced in the following form 



^ D 
B t 

^ V • (p 

0 

II 

T > 

(3J) 


DV =: , 
" Dt 

“ 9 p + X+ 

(7xB)x B , 
0 

(32) 

P 

(9 X 

— r~ 


(X V T) - 9 • (p V ) + rp . 

(33) 


o 
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This sot of equations together with £q« (30) formed a complete description of the 
time dependent problem.. As observed from these equations^ the only electromagnetic 
quantity explicitly involved in this system of equations (30-33) is the mognetic field 
B . 2f B Is given^. this set of equations is reduced to ordinary dynamic equations^ 
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II -3 Disfurbed Thermosphere; Dynom? c Responses of Thermosphere 
II ^3-1 Statement of the Problem 

A sfmpitfied model considers the thermosphere {8Ch200 Km)as a neutral 
inviscid atmosphere with finite electric conductivity * Hail effect can be neglected 
as we have shown in Section II. Thus, the set of fundamental equations is as 
follows: 


+ V ‘ ( pv ) - 0 

04) 

- - 7 p + { j X B ) 

Dt p ^ P 

(35) 

: = div (X V T) + j • E’ 

put 

(36) 

E‘ = E + v X B 

(37) 

V X E = ~ 

(38) 

— * — » — * — ♦ 

j = a(E+ V X B) 

09) 

p = p RT 

(40) 


We further assume that the atmospheric motions are confined to the verticol 
(z ~ axis) and horizontal (x-oxis, £-W) directions. For the magnetic field, since 
observational dataaronotavat!ableat this stage, we assume quite arfaltrofily that 

B (x, y, t) = b (t) exp {- P r ) , (41) 
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where r yx + y and 3 is an arbitrary constant, which Is the measure of 
the damping rat© of geomagnetic variations. This assumed expression for the dis- 
turbed geomagnetic voriotions is such that the maximum disturbance is at the center 
of the region in which we are interested, and it will decay exponentially according to 
the distance from the center, which also implies that the geomagnetic disturbances 
are confined in a finite region. Furthermore, the altitude dependence is also 
Ignored in the present stydyo The quantity b (t) is taken from ground station magnetic 
data. In this calculation the May 1967 geomagnetic storm data recorded at College, 
Alaska; Dallas, Texas; and Honolulu, Hawaii were used. The 1964 Jacchla model 
atmosphere was used to calculate the pre-storm conditions of the atmosphere. The 
set of equations (34-41) was integrated numerically by using a finite difference tech- 
nique. The numerical process is explained in Diagram (2)o It includes the following 
steps « 

1. Calculate E from Eqs.(38)and (41). 

I 

2. Obtain E from Bq. (37). 

3. Calculate] from Eq. (39). 

4. Integrate Bq. (35), determine v for next time step, v (At). 

5. Integrate Bq^ (36), calculate T (At). 

6. Calculate p (At) from Bq. (40), 

7. Determine p (At) from Bq. (34). 

8. Use V (a^) to calculate E (^t) from Bq, (37). 

n-3-2 Results 

Numerical results are presented for winds, joule heating and electric field at 
the center of the storm and an altitude of 140 Km Moy 1967, The calculations 
ore based on the observational data recorded at College, Alaska; Dallas, Texas; 
and Honolulu, Hawaii . The Alaskan results show that the geomagnetic storm can 
generate the horizontal wind (East -West direction) (Figure II -4) on the order of maxi- 
mum 1000 M. sec.“^ and vertical wind (Fig. II-5) approximately an order of magni- 
tude smaller, which agree with the recent electrojet observation Fees [l3j ♦ The 

-3-1 

joule heating (Fig. 11-6) is about a few erg -m -sec. and the electric field (Fig. II-3) 
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on the order of 100 mV * m / which agree well with the results given by Cole [14] 

and Wu/ Matsushita, and De Vries [15]. However, the calculated temperature in the 

present model is rather high. We believe this is due to the fact that the present 

model has not token into account the effects due to viscosity. Considerably smaller 

winds, joule heating, and electric field result when the Texas and the Hawaiian data 

are used instead of the College, Alaska data. The horizontal (Eost-West ) wi nd over 

"1 

Dallas, Texas (Fig. II'IC) is calculated to be ^100 m. sec. , and the vertical 

wind Vx (Fig. 11) is ^10 m. sec. The joule heating (Fig. 11-12) is approximately 

-2 “3 -1 “1 

10 erg - m -sec. ond the electric field (Fig. H-9) is on the order of 10 mV“m 

The Hawaiian results show that v~10m-sec. ^ (Fig. 16), Vy^m-sec. ^ (Fig. 11-17), 

-4-3-1 -1 

Q 10 erg - rti -sec . (Fig • 11-18) and E ~ 1 ni - sec . O^ig . 11-15). The 

results show that the strongest thermal and dynamic responses in the thermosphere due 

to the mognetic storm are in the polar region and they become weaker as latitude 

decreases. 
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n=’4 Enclosure 


In this study, e theoreticol model for the dynamical responses to the geomagnetic 
storm in the thermosphere is derived from magnetohydrodynamic theory. The 
validity of this proposed model is based on the assumption that the thermosphere behaves 
as sn eiectric conducting fluid. This is true because the degree of ionization in the 
thermosphere is 10 ond thus the dynamical properties in the thermosphere are domi- 
nated by the neutral gas. In order to test the proposed model, a numerical example 
is presented. In this caiculation, we have neglected the viscous effects and fixed our 
attention only on o narrow region of the thermosphere, ond constont transport properties 
are assumed In the model. In general, the results we obtained are in good agreement 
with the wind measurements from an Agena satellite [13] and joule heating obtained 
by Cole [3 ] « However, the calculated density and temperoture in this model are 
unreasonably high. We believe this is because we have ignored the viscous effects, 
the altitude dependence, and the gravitational wave effects. Therefore, we shall 
recommend that this model be improved in the following manner: 

1) By including the viscous effect. 

2) By including a global calculation model using a suitoble spherical coordinate 
system. Thus, the latitude dependence can be incorporated Into the model* 

3) Currently, our results are obtained for a point. We hope to do a calcula- 
tion to include whole region of the atmosphere from 90 Km - 250 Km . the 
structure of the winds and electric field con be obtained, 

4) By examining joule heating and the effects on gravitational waves* 

These will form a bosis of our future studies. 
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CHAPTER III 

JU UL 

PROPAGATION OF HYDROMAG NETIC WAVES IN THE IONOSPHERE 
111“ 1 Introduction 

How hydromognetic woves of solar origin heat and accelerate the solar 
wind has been one of the most exciting topics in solar physics in recent years 
(Hundhausen, 1970; Pbrker, 1971). In particular, the dissipation of magneto- 
acoustic waves of solar origin by the various transport phenomena in Coulomb 
collisions may be a major source of heating of the electrons in the outer solar 
corona, within five solar radii, which is a region of transition for electrons 
from a collisionol to a collisionless regime (Hung and Barnes, 1972a; 1972b; 
1972c). A similar procedure may be followed in studying the physics of the 
topside ionosphere. 

In general, propagation of hydromagnetic waves is modified by tronsport 
phenomena due to Coulomb collisions. Collisionol effects can vary from region 
to region in the plasma of interest. In particular, waves of a given period may see 
one region of the plasma os collisionless (in the sense that the wave period is 
short compared with the coulomb collision time) and another region of the plasma 
as collision-dominated. In studying the topside ionosphere, waves of period 
T ~ 10 seconds are of special interest, because 10 seconds is the period of the 
typical ionospheric noise.* Electron collision time at a height of 600 Km is 0.1 
second, and at a height of 1000 Km is one second (Hanson, 1965). The ion 
collision time at 600 Km is 60 seconds (Hanson, 1965). Thus, for wave period of 
10 seconds, in the altitude from 600 Km to 1000 Km, electrons are in a transi- 
tional regime between collisionol and collisionless conditions while ions are In a 
collisionless regime. 


Typical ionospheric noises such as hiss, whistler modes, micropulsations, etc., 
are those which have wave periods of the order of severol seconds, while the ion 
gyrofrequency, under 3000 Km height, has a wave period of several milliseconds 
(Barrington and Fejer, 1965). In the present case, we chose the wave period of 
hydromagnetic wave os 10 seconds because it is believed that those ionospheric 
noises might be relevant to hydromagnetic waves. 

**This chapter was written with Dr. R. J. Hung. 
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It is therefore of considerable interest to investigate how hydronmognetic 
waves behave in this intermediate situation, and in (XJrticular, how they are 
dissipated and what can be the criteria of stability, if instability occurs. In 
this study we consider small amplitude waves in the situation in which electron- 
ion and e! Qctron —electron collisions (but not ion— ion collisions) must be token 
into account. It is a reasonable approximation to assume thot electrons follow 
fluid equotions while the ions ore adequately described by kinetic equations with 
ion-electron collisions. 

In the lower altitude region of the ionosphere, sound waves ploy a significant 
role in modifying the atmospheric motions. Sound waves are important only when 
the wave frequency is much smaller than the collision frequency. As altitude 
increases, the total number density decreases and collision frequency decreases. 
When we reach the point where the collision frequency becomes smaller than the 
wave frequency, we are in the region where the ion--acoo$tic wave, which is the 
extreme cose of the magneto-acoustic mode when the collision frequency tends to 
zero, becomes important. The dissipation mechanism affecting wave propagation 
also changes from transport effects due to Coulomb collisions for the collisionol 
regime to Landau damping for the collisionless regime. 

It is known that three different modes exist in hydromognetic waves, namely, 
Alfven, fast and slow modes (Landau and Lifshitz, 1960). In the laboratory, only 
fast and Alfven modes hove been observed (Schotzmon, 1965). The nonexistence 
of slow mode in nature is because the damping rate of the slow mode is almost 
comparable to the wove frequency, so that in the laboratory there is no way to 
observe the slow mode before it is domped (Barnes, 1966). Thus, in this study, we 
ore interested only in the less damped fast and Alfven modes. 

There is another long-period-wave with a period from 30 minutes to several 
hours called a gravity wove (Hines, 1960). The situation for gravity waves is quite 
different than for hydromognetic waves. The propagation of a gravity wave In the 
lower ionosphere (mostly D and E region) is accompanied by ihe propagation of a 
pressure wave. Therefore, the condition for the existence of gravity waves is the 
same os the condition for the existence of pressure waves; namely, the collision 
frequency is greater than the wave frequency. Thus, in the present study, gravity 

waves are less important. 
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As shown In Sec, II, fhe electrons behave as a fluid whose heat balance 
1$ controlled by thermal conduction and viscous dissipation; and the Ions behave 
much as If they were In a perfectly collisionless plasma, with a slight modification 
due to ion-electron collisions. In Sec. Ill and IV, we derive the dispersion relation 
which characterizes the behavior of hydromagnetlc waves, and several special cases 
of wave l^opagation without using numerical calculation. The criteria for mirror 
and fire-hose instabilities are discussed for anisotropic ions In Sec. V, In Sec, VI, 
we discuss In detail the dompi ng of mogneto-acoustic woves under an assumed 
condition for the upper F2 region; the Ions are heated by Landou damping in much 
the some way as In a perfectly collisionless plasma, and electrons are heated by 
conductive dissipation of waves at about 600 Km altitude and viscous dissipation of 
waves at about 1000 Km altitude. The implications of turbulent or irregular fluc- 
tuations of fields due to mirror and/or fire hose instabilities are also considered in 
Sec. VI. 

The theory of the present paper breaks down when the wave period is shorter 
than the ion collision time or the wave period is longer than the electron collision 
time. The former condition will occur (for 10 second waves) at an altitude lower 
than 500 Km where both electrons and ions are governed by fluid equations. The 
latter condition may obtain In the ionosphere at altitudes higher than 1000 Km 
where the fluctuations in electron and ion components follow the kinetic equations 
for a collisionless plasma. Both of these limits will be analyzed in subsequent 
papers. 


45 



Ill “2 Bessie Equations; A Mixed Model 

As mentioned above, in the upper F2 region, at altitudes of 600 Km 


and higher, electron-electron and electron-ion collisions con affect hydromag- 
netic waves with periods of 10 seconds. We idealize this region for the time 
scale of interest by treating the electrons as a fluid and assuming that the ions are 
adequately described by a kinetic theory that neglects ion-ion, but not ion- 
electron collisions. We further idealize the plasma and magnetic field by 
assuming that they are, on the overage, uniform throughout an effectively infi- 
nite volume, that there h no net average electric current, and that the average 
electron pressure tensor is isotropic, but the ion pressure tensor is not„ Let n, 

T, Pj q f ond tt denote number density, velocity, temperature, pressure, heat 
flux, and viscous stress tensor, respectively, and let subscripts and superscripts 


e and i denote electrons and ions, respectively. Let ^and^ be the electric 
and mogneric fields, e the Ion churge, electrun mosa, und c 

the speed of light. Then, the fluid equati^s for the electrons, and the kinetic 
equation for the ion velocity distribution f. may be written as follows 
(Broginskii, 1965): 
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where P 




dt 




i V 


is the 


electron collision time, and R " - R. is the collisionol momentum transfer 

‘^•1 

from ions to electrons. R is composed of o frictional force R and a thermal force 
R ,R”R+R .The electron thermal flux is composed of analogous 
parts, £ ^ ' Finally, of course, the electromagnetic fields must 

satisfy Maxwell's equations. 

In the present anolysis, we are concerned with waves whose circular fre” 

quency u> is small compared with the ion gyrofrequency, and whose wavelengths 

are long compared with the mean ion Larmor radius. We neglect the frictional 

force R becouse R is important only if t e 1 (condition Ti t 

^ g ~ u e e 

> > 1 is always true for upper F~region ionosphere at 600 Km or higher), where 


' eB/m c is the electron gyrofrequency (cf. Hung and Barnes 1972a; 1972b). 
e e 

The momentum transfer due to collisions can be reduced to 


R = ~ 0.71 n V T 
~e e e 


(2.5) 


where the subscripts II, X, refer to the magnetic field directions ^ / I S 

Similarly to the present order of approximation, the electron heat flux is 


“ K 


V T 

e I ( 


(2.6) 


where K is the coefficient of electron heat conductivity. Furthermore, the stress 
e 

tensor (Braginskti 1965) is 

Here the unit of T is the erg. 
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- e ® 

TT = - ri : w 

fti ft:! «S! 


(2.7) 


where the rate of stroin tensor, w , Is 


(x0 ;^x 


( 2 . 8 ) 


and Ti IS the tensor coefficient of electron viscosity, which is the function of 

ft=S 

a and T and subscripts a and prepresent the coordinates. Under the present 
e ® 

condition Q t the stress tensor has the following form in a coordinate 

e e n V 

system with z-oxis parallel to the magnetic field (Hung and Barnes, 1972b; 1972c): 


n ® 211® — 

zz VO n z 


(2.9) 


(2.10) 


e e 

n ~ rr 
xy yx 


n ® ~ (9^ ( l/Q T ) 
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where the zeroth-order coefficient of electron viscosity is 


Ti® = 0.73 n T T 
VO e e e 


( 2 . 11 ) 


ni-2-1 Electron Dynamics 

In order to solve the fluid“lil<e equations, we assume 

T - T> t 6 t (x , t) 

n ^ n > t 6 n (x , t ) 

p <p> ♦ 6P (x , t) 

B ” <B>+ 0 

E *E (k , t) 

V® 6 v (x , t) 


( 2 . 12 ) 
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where ^ ^ denotes/ ensetrible averaging* We consider the limit of small ampli- 
tudes fluctuations, | 6 / <n > | « 1, etc. If the fluctuations are sinusoidal, 

n • 

i . e./ proportional to exp [i ( k * x - uut -3 where, i = >/^, _&|s. (2.1 - 2.3) 
become after linearization (and neglecting terms of order m^ / m. ) 


6n =<n > (k, 6v ® +k,, 6y ® , 

e e X X I I z 
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Here, without loss of generality we have assumed k (k,, 0, k,,). Combining 

~ X M 

Eqs . (2.13) and (2* 17), we have the following relation for pressure and number 
density 
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To investigate the physical significance of Eq- (2.19), let us examine 
the parameter f' in the following way: 


lim f. 

K -0 r = 4- 

e «j 

(2.2Ia) 

lim 

K “* ® r = 1 

(2.21b) 


e 


This means thot in the two extreme cases of zero and infinite thermal conduction/ 
the electrons fluctuate adiabaticolly and isothermally , respectively. 

The velocity fluctuation of electrons can be easily found in terms of ^ £ 
by usi ng Eqs . (2 . 1 3 ) a nd (2 . 1 8 ) i n Eqs • (2 » H " 2 . 1 6 ): 



where the components of the mobility tensor are 
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is in general a complex funcHon of and a ( 2"^T ^ / m ) is the 

e e 6 

electron thermal speed. From Eqs. (2.23), (2.20) and Broginskii 's (1965) expression 


for K , it follows that 
e 


m. 


(2.24) 


where 0 = 8ti<n><T >/<B>^, u=iu/|k,,|r. , ond C.=<B>/ 

e . e I I A A 

(4 Hm. <n >) Alfven velocity- In deriving Eq- (2.22) we hove exploited the 

fact that I | 1. 

e 


IIH2-2 Ion Dynamics 

Next , we consider the ion kinetic equation, with the object of finding 
an expression analogous to Eq. (2-22), for ions. The first term on the right hand 
side of Eq. (2.4) affects the evolution of the velocity distribution on the time 

1/2 

scale T m, / m t. (m. /m ) which is long compared with time scales 
e I e I I e 

of importance for the wave, and may therefore be neglected. As we pointed out 
earlier, the frictionol force is also negligible. Hence, from Eq. (2.5), the right 
hand side of Eq. (2.4) is just 
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(2.25) 


51 



Linearizing Equation (2.4), we obtain 
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whose solution is 
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Faraday’s law 
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bqs« (2.13-22) permit us to write 6 B and 6 T as linear combinations of 

— e 

the fluctuating electric field components. Then Eq, (2,28) can be solved by 
standard techniques to give an ton mobility tensor analogous to the electron mo- 
bility tensor of Eq, (2»22). However, it is simpler to exploit the fact that Eq. 
(2.28) is equivalent to the linearized Vlasov equation with 6 E replaced by 

6 £ + i (0.7 1 / e) k , , 6 T e t,* Hence, if 6 T is expressed in terms of 6 E / 
^ 1 1 e e ~ 

we moy easily find the ion mobility tensor from the Vlasov mobility tensor. In the 
present case, we have 
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Then, proceeding as In Hung and Bornes (1972o/ b, c), we find 
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where the mobility tensor with anisotropic ions are 
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ijj 


/ ( I k|j I ojj ) , and Z. " Z (y.)and Z.‘ are the pasma dispersion 




function and its first derivative (Fried and Conte 1961). In deriving Eq. (2.31), 
we have neglected the terms proportional to A Z. ' because 
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( 2 . 32 ) 



ni**3 Dispersion RetoHon 

For fhe limit of hydromagnetlc waves, the displacement current 1 
negligible. Then the electric current density may be written as 
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4 - TT 
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where the dielectric tensor 
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From £qs. (2.22), (2.31)ond (3.2) we have 



L 


where the components of the dielectric tensor are 
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The requirement that Eq- (3.1) be consistent with Maxwell's equations gives 
the dispersion relation 



Substituting Eq* (3.3) in Eq. (3.4) gives 
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to lowest order in uj/q. . The first factor is the usual Alfven wave dispersion 
relation 



(3.6) 


and the second factor, which gives the magneto ’•acoustic dispersion relation. Is 
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By using Eq. (3.3), we can rearrange Eq. (3.7) by straightforward calculotion 


to give 


[ ] + ] cot^ 0 - S (u, ) (3.8) 


where 
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and u = (U / ( |k j| I ) 


(3.15) 


We notice that if dispersion relation (3-8) is satisfied for a wave ('i>, k ) , 

it is also satisfied for a wave propagating in the opposite direction, I.e., for the 

woves { ^ " Ik ) ond ( “ ^* # ) • This is so first, because (3.8) is invariant to 

the transformation ) and secondly, because the transformation 

(oj, Ik) - (-UJ* , k ) causes , C -- C* , T T * , Z' Z'*, 

2 2"^ IClLv vee li 

u ( ^*) / so that Eq . (3.8) is transformed into complex conjugate. Hence, in 

particular, we may always choose Re r\ > Q and ReC > 0 . 

TC ~ V ~ 

Note the superscript * indicates the complex conjugate. 
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in -4 special Coses of the Dispersion Retotlon 

Eq- (3-6)9ives tl^e dispersion relation for Alfven waves in a plasma with 
anisotropic ions* The analog of Eq. (3.6) in a purely collisionless plasma is 






( ^ - 3 „ ) 


(4.1) 


where 


“ r- / J 1 \ ' 


\ 6 ^ 
. For isotropic electrons, < 


p® ~ 0^ , Eq. (4.1) is reduced to Eq. (3.6)* If the ion collision frequency is high 


enough to make the ion pressure tensor isotropic, the dispersion relation (3.6) becomes 
identical to the transverse Alfven wave of magnetohydrodynamics. 

The mag neto~acou Stic dispersion relation (3.8), which is much more com-* 

plicated than either the coilisionless or collision-dominated cases, is transcendental, 

and must generally be solved numerically. However, we first consider a few in” 

teresting speciol limits that can provide insight into the character of the waves 

without numerical calculations . In the present case, S depends, in general, on 

'JDand 0 through and C which are functions of the electron-electron 

. TC v 

collision time T , and T depends on the altitude of the ionosphere. 


in~4-1 Propagation Parallel to 
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Eq, (4»2) shows fhat one of fhe mognefo-acousHc modes propagates as an Alfven 

wave when 9 = 0 in a plasma with anisotropic ions • Eq. (4.3) corresponds to a 

compressive wave propagating along < B ^ « The analogous wave in collisioncl 

magnetohydrodynamics is the sound wove# and in a coUisionless plasma the analogous 

wave is the ion-acoustic wave, which generally undergoes Landau damping (Fried 

and Gould, 1961). When ^ is on the Order of or smaller than unity, which corresponds 

e 

to the altitude of about 1000 Km in the ionosphere, C is also on the order of or smaller 

v 

than unity, so It loads to 


I T] I > > 1 and r 1 . 
TC e 


(4.4) 


This means that when the electron collision frequency is on the order of the wave fre- 
quency, or at dn altitude of about 1000 Km, the viscous dissipation is more effective than 
conduct; ve dissipation . under this condition, tq* (4*3) becomes 



2 < T .1 > 

, . 4 ^ 

1 + i C 


Z‘ 

s 


3 V 

(4.5) 

<T > 
e 

1 • 1 4~ ^ 1^ 

O V 

If Eq. (4.5) has any weakSy damped solution, 

F i.e., solutions such that 

1 Im'U / 

Re 1 1, and S S 

V 

0 

1 - R e S , we 

V V 

would have 


ReZ.' ' 
1 

2 <T ' > 

II 

1 

(4.6a) 

< T > 

e 


Im Z.‘ 
1 

2 <T ‘ > 

11 

4 ^ o 

3 ' V 

(4.6b) 

< T > 
e 

1 ‘ (4- ^ 

O V 

Since I m Z.' 5 0 for rea 
1 1 

1 arguments, and Eq. (4.6b) cannot be satisfied for 

1 I m u / u I < < 1 unless 

Re u > > p 

i 

II 

(4.7) 
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Therefore, Eqs. (4,6a) and (4.6b) become 



Im UD 
ReUJ 


0.49 T 

6 


which ore valid only if 

3 >> 3 . 

e I 


(4.8a) 


(4.8b) 


(4.9) 


This condition is nothing but < T > >><T. > which is just the condition for weak 

e « 

damping of ion acoustic wave in o perfectly collisionless plasma (Fried and Gould, 

1961). 

Next, let us consider the other limit, i .e-, when is on the order of or 

smaller than ( m / m. ) , which corresponds to an altitude of about 600 Km or 
e » 

higher in the ionosphere- Under this limit, conductive dissipation becomes signifi~ 
cant. Therefore, Eq. (4.3) becomes 


Z.’ 



< j > 

e 


15 + 4 I 

25 I 4 I 



(4.10) 


where “k,, K /u;<n >. For the weakly damped solutions, i.e-, 

TC II e e 

i I 1/ and t) p=i ~ ^ ^ TC ' would have 



15 H 4 1 I ^ 
25 » 4 I 


(4.11a) 
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^ ‘it 

<T > 
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4 X] 
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25+ 4 h 


(4.11b) 
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Using the same argument as in Eq. (4,7), Eqs. (4.11a) and (4.11b) become 


(4.12a) 


(4.12b) 


which are valid only if 


2 11 ° 
'tc 

15 ' 

< < 1 

(4,13) 

Parti culorly^ in the cdiabctic limit 
in the ionosphere/ 

r|° -*0* corresponding to t- 

1 

-region or lower 


<T > 
e 

(4.14a) 

( h j - 3 

t 

m. 

\ 

Im'JJ 

2 _o 

(4.14b) 

Re^ - 

15 ^TC 

ore sonic waves propagating parallel to < B > . In the isothermal 1 
-♦ «> / it corresponds to 2000 Km altitude in the ionosphere, 

Tmit 

/ . 

A 

fv 

(4.15a) 

\ ) 

- m. , 

1 

Im 

~ - -1_ ri° 

(4.15b) 

Re 

2 ^TC - 


which resuifs in an equation similar to Bq, (4,8) although the damping rates are 
di fferent . 
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In summery, we have shown that viscous damping is significant at the higher 
altitude (about lODO Km), and conductive dissipation is overwhelming in lower al" 
titudes (600 Km or higher) of the ionosphere. In practice, oil ion^coustic waves 

propagating parallel fo ^ above the F region, and sonic waves at E-region or 
lower, are dissipated, only if that the case of < > >> <T. > will assure the 

propagation of the weakly damped ion-ocoustic waves. 


ni~4~2 Propagation Transverse to <^> 

I 

When kX'^B^, Z “*0 and 0 . Under this condition, i' 

~ ~ I TC e 

becomes 5/3 and expression (3.8) simplifies considerably. 


tel 


^ i A 

X e \ 6 


(- 


■) 


(4.16) 


Therefore, we have 




2<T > ^ ^ 

c 2 , . 1 

A m. 3 m. 


C <T > 

V e 


3 m. 


(4.17) 


If C is small , and C “ReC ^^ImC , then the real part of (4 . 17) is 

V V V V 




2 <t' > I i <T > 

c2 , i_Lt — i_X_ 

A m. 


(4.18) 


and 


I m 

Re ~ 


V k 


<T > 

e 


^ (R e uu ) ^ ""i 


(4.19) 


Eq. (4. 18) is intermediate between the analogous expression which arises when ions 
are colli si onal 


(-f) 


2 ^ . < T. > 1 <T > 

„ ^ 2 o 1 e 

A 3 m. 


(4.20) 
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and that which arises when electrons are collisionless 




<T^> 


+ <T, 


m. 


(4.21) 


In a coil SSI on! ess plasnrta, strong electron Landau damping can occur for 8 slightly 
different from tt/2 (Barnes 1966). However, when k Is exactly perpendicular 
to<B >, this Landau damping vanishes. At the lower altitude of the F-region where 
conductive damping overwhelms viscous damping, conductive dissipation vanishes 
at 0 tt/2. In the case of a higher altitude, say 1000 Km, where the viscous 

damping becomes significant, the result is different. No matter how small the co“ 
efficient of electron viscosity Is, as long as is not zero, viscous dissipation 

always can be found even when k is exactly perpendicular to ^ • 
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Ilia’S Criteria for Mirror and Fire-Hose Inst-gbilUies 

Most of the known plasma instabilities, which are termed oversta bill ties, are a 
kind in which iw is complex rather than pure imaginary (Stix, 1962)- The non- 
occurrence of overstoble solutions has been shown by Vedenov and Sogdeev (1958)/ 
Chandrasekhar et al (1958), and Barnes (1966) for the case of purely collisionless 
plasma. The significance of the nonexistence of overstable solutions is simply that the 
condition of marginal stability for the plasma is given by the dispersion relation 
with " 0. This conclusion remains valid for the case of collisional electrons and 

collision! ess ions. 

Let us rewrite dispersion relation (3.8) 

D (u, 8 , ^ ^) “ 1 < ^ ” ^^11 ^ ~ cos^ 9 - S (y>, t^) sin^B 

L J ® (5^1) 

Assume that and 0 , or u ond 0 are an overstable solution of Eq. (5.1). Then 

V 



(5-2) 

Now 

Sgn 1 Im (Zj ) ] - Sgn (u ) , 

Sgn 1 .! Sgn I ^ Sgn ( u) 

where 

Sgn 1 X I X / I X I , 
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so 


that 


Sgn [ Im D ( u, 9, t ) ] “ sin ^ j ^ 


m! 


II 


Im(z;) 


/ 61 ' 

C ll-s- 1 

V \ Ct 




’ r 


8 r e I + 4 g r ^ 

e X e e 



(u) 


(5.3) 


Hence/ if sin 9 ^ 0, Im (D) cannot vanish unless u - 0, so that Eq- (3«8) has 
no overstable solution unless sin 9 ~ 0- If sin 0 0, the dispersion relation re- 

duces to Eq • (3 .6) . 

Tho nonexistence of overstable solutions recovers the familiar fire-hose and 
mirror instabilities for collisional electrons (isotropic pressure tensor) and collision- 
less ions (anisotropic pressure tensor). The criteria of the instabilities in the present 
case can be expressed as follows: 


1 + 



(5.4) 

(5.5) 


Eq- (5.4) is the Alfven fire-hose instability to which Eq. (3.6) is relevant, and Eq. (5,5) 

is the magneto-acoustic fire-hose Instability to which Eq. (3.8) is relevant. In both 

coses, there is an angle 9 at which the stability is nrKarginal . For the magneto- 

o 

acoustic fire-hose instabi lity, unstable waves propagate at angles 9 such that 


0 < 9 < 0 or n - 0 < G < n , (5-6) 

— o o “ 

and for the mirror instobi I i ty, unstable waves propagate at angles 0 such that 
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TH© equqtton for 


1 • 


This angle 6 is 


e < 6< JL or — <9 < n - 0 (5.7) 

o 2 2 o 


0 Is just D (0, 0 , 0 ) " 0 , which may be rewritten as 
o o 



(5.8) 


real only if one of the instability criteria (5-4) and (5«5) holds. 
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Ill “6 Numerical Solutions and Discussion 

m mi < I III ■ I , III tm .j m i ■■ ■■ M 

As has been shown in Section III, the dispersion relation which governs the 
behavior of magr^eto -acoustic waves is quite complicated. Unfortunately, analytic 
expressions for the damping rates are tractable only in very special limits. There- 
fore, we have resorted to numerical calculations to compare the damping rates of 
hydromogneti c waves for various values of the collision frequencies which character- 
ize the altitude of the ionosphere. 

We are interested in values of electron collision frequency appropriate to 

the ionosphere, at a height of 600 to 1000 Km. As jointed out earlier, that transition 

periods of electron collisions, Ui t , from the order of m /m. to unity or less 

approximately correspond to the altitudes of the ionosphere for the present study. 

Furthermore, the present theory is of special interest because it clarifies both the 

charocteri sties of the upper F2 region which extends to the magneto-sphere and the 

transition between collisional and colilsionless behavior of . hydro magnetic waves. 

Therefore, we consider solutions of the dispersion relation over a range of uu t , 

0 

ranging from UU T m /m. up to uu T '>^0 (1) for low 3 (r:( 0. 4) pi asmo. 
e e 1 e 

For the isotropic pSosma ( < Tj|> < Tj^ >), o), the numerical results 

of the dissipation rate as a function of the propagation angle 6 are obtained for 

higher ( uj T “ 0.083) and lower ( ojt ^ 0,46), electron collision frequencies 
e e 

at 8 ^ 0,4. It shows (See Figure) the present results are similar to those obtained 

by Barnes (1966) for the collisionless plasma, therefore, we may conclude that the 

two peaks appearing in our calculation are due to ion Landau damping for 

0 50°, and electron thermal conduction and viscous effects for 0 > 50°. 

Hence, the electron heating is always more important than the ion heating under 

the present conditions. This result explains the facts that the dissipation of 

hydromagnetic woves might play an important role in creating the phenomena 

< T > > > < T. >, in the upper F2 region, which agrees with observations 
e I 

(Comfort, 1970). 
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For the plasma with anisotropic ions, < T jj > |^ < Tj^ ># the mirror and/or 
fire-hose instabilities might occur bosed on the instability criteria (5.4) and (5«5). 
Although the experimental data are still not enough to prove how far the anisotropic 
ions can sotisfy the instability criteria, it is expected, however, that mirror an6/or 
fire-hose instabilities might be relevant to irregular fluctuotions or turbulent fields 
observed in the upper F2 region. 

In summary, the present theory reveals the following properties of fast-mode 
magneto-acoustic waves under conditions expected in the upper F2 region: (1) 

Conductive domping is dominant at about 600 Km altitude and viscous damping is 
overwhelming at 1000 Km altitude; (2) Dissipation of hydromagnetic waves creates 
the bi-thermal phenomena (< T > >> < T. > ) which assures the propagation of 
ion-acoustic woves in the ionosphere; (3) The irregular or turbulent fields observed 
in the F2 region might be caused by the instabilities due to anisotropic ions (mirror 
and/or fire-hose instabilities); (4) Conductive and/or viscous dissipation is more 
effective than Londau damping due to wave-particle resonance interaction in the 
ionosphere, while Landau damping is the primary dissipation mechanism in the 
interplanetary space. 
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Fig. III-l Damping rote of magneto-acoustic mode as a function of propagation 
direction for higher (tu = 0.083) and lower (u) “ 0. 146) values of electron 

collision frequencies in an isotropic plasma ® ^ ~ 
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CHAPTER IV 

INTERPLANETARY SPACE ENVIRONMENT 


IV- 1 InlToducflon 

It is ’well known that ordinary gasdynamic theory has heen used 
successfully for a zeroth order theoretical continuum approach to 
the prohlem of solar wind interaction with Earth, and very possibly, 
Venus and Mars (Spreiter and Alksne, 1970; and Dryer, 1970) » The 
interplanetary magnetic field was considered only in an implicit 
manner in the early solutions referenced in these papers. Explicit 
consideration Ox the field - recently completed by Shen (1972) and 
Hirsh and Re shot ko (1971) - shows that the Lorentz force is effective 
only in the neighborhood of the subsolar region. The zeroth order 
approach has also been suggested (in varying degrees of quantitative 
analysis) for the case of moon by Lyon et al. (1967)} Michel (I968), 
Wolf (1968), Siscoe et al. (I989) and Spreiter et al. (1970). The 
boundary condition on the sunlit hemisphere for these continuum 
studies was modified to permit absorption of all incident particles. 

An alternative approach for Moon has been used successfully by Whang 
(1968, 19^9> 1970) who used a ‘'guiding-center” model which reproduced 
experimental magnetic “anomalies” in the woke region (Ness et al. 

1967; and Ness, 1972). This model (also with the complete absorption 
assumption) uses the microscopic kinetic equation to describe the 
motion of guiding centers along, but continuum flow perpendicular to, 
field lines. Whang anti Ness (I97O) further demonstrated the first 
experimental observation of standing magneto- acoustic waves which are 
inclined from the lunar limb in the direction of the "Mach cone." 

In order to interpret the data on a theoretical basis, they, returned • 
to the concept of the anisotropic propagation of magnetoacoustic waves 
as derived not only from the solution of the linearized continuum 
equations but also from the solution (Wiiang, 19?0) of a guiding center 
plasma. This approach was formalized by Wang (1971) who used the 
same absorption boundary condition on the sunlit side of moon and 
found similar density profiles in the wake with a three-dimensional 
kinetic description. In effect, then, moon might almost be considered 
to be no more than a point disturbance in the solar wind. Yet, the 
realities of magnet5.c and plasma density anomalies discussed earlier 
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strongly suggest that moon and other ’'small" obstacles require fur- 
ther studies which start from first principles. For example, one 
might conceivably consider a boundary condition whereby particle 
reflection may be more specular-like, as in the case of an iono- 
pause or magnetopause. 


It has been suggested by Dryer (1970) that the obstacle's 
"Knudsen number , (where r^ is the ambient proton's thermal 

gyroradlus and d is the obstacle’s physical or effective diameter) 
be used as a rough guide in choosing a basic approach to the inter- 
action problem. That is, if r^/d for an obstacle (a comet, say, at 
large A.U.) is large, say >10^, we would likely use kinetic theory 
for the study of its interaction. Conversely, if r^/d were <10*”^, 
we could proceed with great confidence with continuum theory. A 
"gi’ay" area lies somewhere between these somewhat arbitrarily-chosen 
numbers. By analogy with hypersonic flight through the whole gamut 
of Earth-bound atmosjjheric and ionospheric levels, this "gray" area 
could be labeled as a "transition" region which conveniently covers 
our Ignorance of the correct procedure for analysis. Alternatively, 
one might suggest the electron inertial length as the characteristic 
Icngtli of the ambient solar wind. The basis for this suggestion 
appears to be in its I'elation to shock wave thicknesses observed, in 
the laboratory and in space. Another' candidate is the proton gyro- 
radius, based on the directed velocity, which is the characteristic 
thickness of the Earth’s magnetopause. A related length scale could 
be the geometric mean of the proton and electron gyroradii. Detailed 
studies would require consideration of the boundary’s nature: 
magnetic, nonmagnetic or absorbing. 


We could, if desired, extend this approach beyond planets, 
comets, and asteroids (Creenstadt, 1971 ) to include our heliospheric 
motion through the interstellai' iftediiim. It was in the spirit of 
reducing the degree of random speculation to a more systematic theo- 
retical. approach that this idea was proposed. It is the purpose of 
this paper to continue in tMs spirit by providing an alternative 
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approach to our "shopping list" of theoretical approaches. Thus, 
we extend the well-known kinetic theory for particle interaction with 
a sphere (Karejncheti and Sentman, 19^5; and Sentman and Karamcheti, 
1969) to Include the effect of an isotropic and anisotropic potential 
to represent the surfaces of some obstacles . The result may then be 
diffuse or even specular reflection. We are still interested, 
regardless of tne outcome of the point regarding the boundary condi- 
tion (on the sunlit side, say, of the moon) in the limiting case of 
perfect absorption (c.f., Wang, 197I) . This can easily be specified 
by letting reflection be eq.ual to zero. We will include it, however, 
In order to maintain generality. Thus the applicability of the 
boundary condition assumed herein to a specific obstacle (say, moon . 
Mercury, etc.) is not our immediate objective. It is our purpose, to 
demonstrate a technique which, in principle, can be adjusted to in- 
clude other boundary conditions and the interplanetary magnetic field. 
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IV ■'2 Discussions of Observations 

Excluding ionospheric satellites from our discussioUj the only 
potentially ’'small’' obstacle for which observations exist is the 
Moon. Some recent experimental observations on and near moon may be 
relevant within the context of the fundamental kinetic approach which 
will be discussed in the next section. Other observations with 
respect to a possible limb shock (or Mach wave) as .well as the mag- 
netic wake anomalies have been extensively discussed in the litera- 
ture (c.f., Ness, 1972; and Sonett and Mihalov, 1972) and, therefore, 
are not repeated here. 

Figure 1 (based on Fig. 1 from Freeman, 1972) schematically 
shows the Moon's position with respect to a nominal position for 
Earth’s shock wave on several dates;' 13-15 December I 969 and 
27 March 1970. On the latter date, Neugebauer et al, (1972) reported 
solar w'ind observations obtained with a proton spectrometer which 
was part of the ALSEP 12 instrument package on the lunar surface. 

As noted in the figure, the package was on the sunlit side of the 
Moon (with instrument normal directed about 35° to the east of the 
sun) and in front of the Earth's shock on 27 March 1970* Figure 2 
(from Neugebauer et al., 1972) shows the velocity and density of 
the solar wind as measured on the surface of irtoom^ 

Vela 5 find Ogo' 5 were also in front of the bow shock; the 
approximate position of the latter spacecraft is shown in Figure 1. 
Their measurements are also shown in Figure 2 In order to present 
"freestream” conditions. The time prior to 07 II UT is of interest 
here. (An interplanetary shock was measured at that time. The 
data subsequent to Oyil UT are of secondary interest here except for 
noting the related shock jump conditions at all three instruments. 

The time scales for these two spacecraft were shifted so as to coin- 
cide with ALSEP 12 at the shock wa,ve.) Note that both spacecraft - 
in the aj:ribient solar wind - measured velocities and densities which 
were, respectively, higher and lower than those measured on the lunar 
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surface. Neugebauer et al. (1972) suggest that this compression 
effect vas caused by local lunar magnetic fields (c.f., Sonett and 
Mihalov, 1972). It is suggested, however, that the results discussed 
later in terms of density contours present a physic ally -plausible 
alternative explanation. 

Several other experimental observations associated with the 
moon may also be relevant to the present study. Figure 1 also shows 
the Moon's approximate position on 13-15 December I969 when the 
ALSEP 12 site was deep in the lunar night. Freeman (1972) has re- 
ported persistent observations of positive ion bursts, "probably of 
solar wind originlV, when his suprathermal ion detectors at ALSEP 12 
and 1^ are located on the night side as indicated for the particular 
case of ALSEP 12 in Figure 1. The low flux suggested later by den- 
sity contours and velocity vectors behind a typical spherical object 
may provide an explanation for these observations , 

Another set of particle observations in front of the moon has 
been reported by Prakash (l970) . Figure 3 shows contours of relative 
2-5 keV electron flux on the sunlit side of theo;niooii as measured by 
the MIT plasma cup on Explorer 35 during — 700 orbits (outside of 
the Earth's shockwave). These electrons, of course, represent only 
the high energy tail of the electron distribution function. Because 
of charge neutrality which exists in the solar wind, the data shown 
in Figure 3 must represent only the "tip of the iceberg" of an 
equally dense presence of protons in or near the same locations in 
front of the moon The relationship of these preliminary observa- 
tions to the density contours noted later suggests again that the 
kinetic approach may be physically plausible for the moon. If this 
suggestion is true, the presently- accepted concept of complete solar 
wind absorption on the sunlit side, with little or no "sputtering", 
may require reexamiaatlon.s< 
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IV-3 Analysis 

Let us consider a collisionless neutral plasma (i.e., solar 
wind) which flov/s past any "small” celestial object. On the surface 
of the object is an assumed angularly distributed positive electro- 
static potential field, A schematic diagram of the problem is de- 
picted in Fig. This problem is equivalent to that of a free 
particle flow which interacts with a solid sphere with a body force 
due to the electric field on the surface. Hence, the mathematical 
model for the present problem can be described by the colli sicr.lecs 
Boltzmann equatioii with the body acting as an external force or. the 
gas, thereby eliminating the need for gas- solid surface boundary 
conditions. For simplicity, we shall consider steady-state icr. 
motion w;i^hin the solar wind; thus the governing equation for the 
present problem is 


V ® 




I'l ' 


where f*=f(r, v) represents the singlet local velocity distribut icr, 
function of the ions; v, the' ion velocity; F, the external force 
given by 

£= -grad$(J (2; 


where $(r), a scalar function of r and 0, is the potential fu.ncticr. 
due to the properties of the surface of the obstacle under consider a 
tion; m, the mass of an ion; and, finally, r is the position veer or. 

The solution of the collisionless Boltzmann Equation has leer, 
discussed in detail by Karamcheti and Sentman (I965) who direct 
attention to the advantages that result when solid surfaces ere in- 
cluded- in the kinetic problem as external forces. Thus, the velccir 

distribution function f(r, v) to satisfy Eq. (l) for |r|-« “ car. ce 

* 

written in the following form: 

^ i ~ A exp [ - B (v - U^)® ) 
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with . 


A = n^{m/277kl^)V^ 

B=m/2kT^ . (3a) 


and where r a.nd v are expressed in a spherical coordinate system. 

represents the incoming solar wind velocity, and T^, the ajnbient 
ion temperature . 


The velocity distribution f (r, jv) is calculated numerically 
for each step. Initially, we assign values of the distribution 
function to the grid points in r-v space. We convect the values of 
f according to Newton's law; 


df 

3T 


e 0 


dr 

dT 


r= V 


( 4 ) 


m 


« - grod $ (r) 


J 


Thus, the distribution function can be calculated successively. It 
is convenient to nondimensionalize all the variables by measuring 
all the lengths, velocities and times in terms of the object’s rail 
a; the thermal speed of the particle, c = i/b^; and the time, aB^, 
necessary to travel one ”a” at the particle’s thermal speed, 
respectively. Thus, we ceilculate the distribution function follow! 
the path given by Eq. (h). Tlius , a non- spherically symmetric poten 
tial which represents the body surface can adequately be taken into 
account. It is important to note that the present solution is the 
lowest-order approximation to the solution of the Boltzmann equatic 
and consequently, represents a process which ignores inter-particle 
collisions. Some additional details are given in the Appendix. 
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On the surface of the body, the velocity distribution function 
for the particles emitted from the surface is expressed by 



vhere B = — ), 0 and v are the constants of integration which 

V 

will be determi.ned from the trajectory equation. In this study, the 
accommodation coefficient is assumed to be unity. More details 
concerning ■this point can be found from Prager and Rasmussen (19^7)* 

Knowledge of the singlet local velocity distribution enables 
the calculation of the local density and velocity by taking moments 
with respect to the local velocity distribution. 

1V“4 Method of Solution 

Now, the task left to us is to calculate the velocity distribution 
using Eq, (4). A detailed discussion of this method is given by Prager 
and Rasmussen (19^7) < shall not repeat it here; however, for 
completeness, we shall outline the method briefly as follows. 

We prescribe a boundary condition for the distribution function 
at t = t^ in terms of some given function f^ , then, 

X (5)- 

Thus, the general solution can be expressed in terms of the 
function f^ by 


f(r(t), V (0) = h (l (L' X.) x) '> 


( 6 ) 
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where 



(T) 


which were ohtained by integrating Eq. (^) . 

We «iay note that as t the initial value of f is 

recovered. The vectors v and r represent the six constants of 

o o 

integration required for the equation of motion. Therefore, the 
general solution is an arbitrary function of the constants of 
integration of the characteristics of Eq. However, this arcitrar 

function is selected to agree with the boundary condition. 

For the present problem, the spherical polar coordinate syste:i 
for the physical space is chosei^ then the solution to the collision- 
less Boltzmann equation corresponding to the condition Eq. ( 3 ) is 
given by 

f (r , V ) = A exp (r, v)+ U® 




I 


Now the task is to determine v (r,v), v (r,v), cos 

0 (r,v) and v^(r,v) by solving the equation of motion of a particle 
a a ^ 

in a potential field $(r,0). 
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After some mathematical manipulation, the distribution function 
for the particles coming from infinity and those emitted from the 
surface are now 

f » X) = [ - S ( v= + ^ ’ 2 |v«+ ^ • cos 0 )' 

end 



n ( e ) 

Q 


m 

Trnrr 


Z/2 


exp 


w 


B (v® + il 

w m 



( 10 ) 


The functions cos 0 and Q again will be determined from the 

00 a 

trajectory equation of particles » 

For convenience, we shall write the total local density in the 

form 


n (r, 0) n (r, e ) n • (r, 0) 

J emit 


(11) 


where (r, &) is the net contribution to the number density by the 

incident particles, i.e., the contribution of those particles that 

are moving toward the object with the effect of those particles passing 

by the object; and is the contribution to the number 

density by the particles emitted from the surface with the temperature 

of the body surface, T . This is the effect of the boundary condition 

w 

on the surface. 
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After ao]T\e mathematical manipulation, ve obtain 


O CO CO 


n, (r .0) = 






f, d^^ dXg dX^, , 


( 12 ) 


with x®= v® + — 4 (r, 6) , 

m 

vhere X ^ (X^ / Xq# X^) # v ~(v^ / Vq / v^) 


and 



Kt)’ 


+ 1 ”2 X cos S+ 2 


XgSin 


0 


(13) 



, the particle speed ratio. 


n (r, 9) = 
emit 



dX^dXg dx^ (Us) 

Where n(8^) 'will be determined by the assurance of continuity in the 
characteristic equation where n(9 ) will be found in terms of n , 

a 05 

The mean velocity field is detennined by the first moment of 
the velocity distribution function, 
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Q _ Co 00 


u. (r, 0) = —5 — _ 

> ' ' n(r, e; 





0/2 






,0O CO oo 


*/■£/. ’‘i('*-^-^) 


1/2 


^emit ^ 


J 


with 


j « r, B 


(15) 


Klns-lly^j ths pctcntio-l on "tho surfs,c6 is dsnlvod froin fi. cim!rgGd 
particle interacting with an angularly-distributed charged particle 
layer on the surface and is given as follows; 


I (tf 0) 



Q (cos 0')V^ d 9' 

- 2 (L) cos (9 - 0-)jjy^ 


(16) 


on^/or 

$(,r} ^ ; n = 2 and 5 in the present study^ (17) 

whore Q Indicates a constant representing the product of quantum 
efficiency and subsolar potential and Q’ Is the usual force constant, 
A repulsive electric field of 1 Vm'^( volt/m) at the subsolar point 
will he assumed for both asymmetric and spherically syinmetrlc cases. 
In the latter case, of course, § is independent of 6 so that the two • 
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cases are normalized at the subsolar point. 

In the figures which are discussed below, the nondimens ional 

radial distance, r/a, is denoted as R. Figure 5 shows the result 

of the integration of Eq. (l6). It is seen that ^ becomes essentially 

independent of 9 within several radii, thus, the main effect (of the 

potential which represents the force on the fluid) is confined to 

small values of R where a shock wave could eventually develop. In 

such a calculation, we have weighted the potential energy by the 

reflected particle's thermal energy, kT , as seen, for example, in 

w 

Eq. (10). 

IV-5 Numerical Results 

A steady solar wind is assumed to be flowiiog past a "small" 

object with the following assumed parameters; U = hOO km sec"^. 

“3 * 5 o 

= 5 cm , and T =10 K. The particle speed ratio, S, is 

therefore 9.85. Other typical solar wind parameters will generally 

give S > 1. Sentman and Karamcheti (1959) show that the results 

for density contours are insensitive to variations in S once S > 1. 

We will assume that the present case, with the electric field, is 

similar and will not investigate this point any further in the 

belief that the present case will be sufficiently representative. 

Using Equations (ll) through (l4) together with (if) through (lO), 

the density distribution along the sun-obstacle axis (i.e., 6 = 0°) 

is first found. Figure 6 shows n/n plotted versus the non-dimensional 

00 

radius, R, in the absence of a field on the surface as well as with the 
presence of a field as discussed earlier. A sharp rise from the 
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ambient value of density to 2 and 3 times tliat value occurs within 
half of the obstacle’s radial distance from the surface for the two 
cases ^ respectively. The results for the major portion of the 
affected flow field is shown by contours of constant density (with 
effective non-central force field) in Figure It is seen that the 
surface density rapidly decreases as 0 moves toward tiie object’s 
limb. A substantial paucity of particles is also seen in the 
immediate base region as shown also by Wang (l97l) for the 
alternative case of absorption on the sunlit surface. Of partic-jlar 
interest, however, is the steep density gradient in the immediate 
vicinity of the limb. This result is, not unexpectedly, similar to 
that shown by Sentman and Karmacheti (1969) for the case of no 

The density contours for S « 10 (with and without a 
potential field) are also given by Dryer (1970). Also, the conto'jrs 
of constant density with effective central force field which are 
described by Eq. (17) are shown in Figs, 6 and 9 for n - 2 and 3, 
respectively. Note that there is no substantive difference betv^een 
the case of an asymmetrical force field (which represents a surface 
electric field as suggested for Moon by W80,bridge (1968, 1972) and 
Grobmon and Blank (l9£>9) Q- spherical Inverse force field. 

Finally, Equations 16 and 17 provide additional information ate 
the flow field, namely, the streamline pattern, Figiare 10 shows 
the results of the calculation for non-central force field where 
the length of the individual velocity vectors is proportional to 
tlielr magnitude. Note, for reference, the length of the undisturtei 
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velocity vector at the left side of the figure. The expected 
deceleration of the flow on the sunlit side of the obstacle with 
subsequent acceleration in- the expansion phase around the limb 
and in the wake is clearly seen. Mention has already been F:ade 
above of the observation by Neugebauer et al. (l972) which my be 
relevant to the present study, it is seen that the flow^ on the 
average^ is similar to that for continuum flow around a sphere. 

As noted earlier, the interplanetary magnetic field has been 
explicitly neglected, Fvirther studies (with the present boundary 
condition) stiould, of course, take B into account, but it is possibl 
to make some qualitative observations regarding its effect on the 
flow. When v is not parallel to B, the induced electric field in 
the object’s frame of reference is - v x B. Thus, the Lorentz fores 
is proportional to (v x x B. We refer here to the 'volume current 

«*>u< 

considered in MHD analyses and do not consider charge separation; 
thus, the E X B drift can be neglected. Assuming, then, that 
B = (B^, 0, B^) in a spherical coordinate system whose axis is at 
the center of the sun, and ^ = (u,o, O)^ this force has the follovir^ 
components: (- u B^) f + (u B^ ) cp, Fqr the chosen luidlsturbed 

velocity components and "toward the sun" polarity of the field, the 
lorentz force acts within the ecliptic plane to cause a distortion 
of the flow field. The density contours and velocity vector field 
(Figs. 7 through lO) would no longer be perfectly symmetrical. V.-ar^ 
( 1971 ) has shown, for example, that the density distribut'ion shows 
"complicated periodic structures" for cases of large gyroradli. 
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The direction of the force, indicated above for the example, vill 
favor a higher density concentration on the sunward, morning quedr 
of the flcfw field o The results presented here, then, must be 
considered to be a qualitative estimate of a self-consistent kinet 
calculation which explicitly incorporates the effect of an arbitra 
oriented, frozen-in, magnetic field c The skewing of the density c 
(as well as other parameters) has been explicitly found in the 
continuum case by Shen (1972) and Hirsh and Reshotko (l97l). Sons 
of the latters’ work is discussed by Dryer (1970) . An additional 
manifestation of this effect would be the unsymmetrical "Mach cone 
which has been observed by Whang and Ness (l970) in. the case 
of moon and which is anticipated by the limb density gradient in 
Pigs. 7> a^d 9 for the zeroth order solution. 

3V"6 Enclosure 

Sentman and Karamcheti (l9&9) have pointed out that the spaci 
of the constant density contours in front of the sphere Is relativ 
insensitive to varlattons in the speed ratio. Sconce S > 1. Behind 
the sphere, the contours (hence the wake) elongate while there is 
a decrease of the width of the compression . and the "expansion fan” 
which extend rearward from the limb. Althoiagh results (which incl 
effects of the electrostatic potential) for S different froin 10 ar 
not presented here, we believe that they would produce similar eff 
as noted above. It is seen that the ridge of increased density 
^ ^a>) - fo reflected protons - extends downstream In the 
direction of what might be called a developing ‘’Mach coney The 
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results suggest, for high Knudsen numbers, that effective collislonal 
encounters would first become important along this ridge of high 
density gradient.. Sentman and Karamchetl. (I9b9) have suggested 
that" perhaps it is the high-density ridge that develops into 
a shock wave as collisions become important" Also, the shock may 
indeed form in a high Knudsen flow, but it would probably first 
appear dfjwnstream from the body and, as Sentman and Karamcheti 
point out: " ... the stagnation streamline would probably be the 
last region in the flow to become aware of the presence of the 
shock"! Flybys of objects (such as Mercury and Pluto which are 
potential candidates for high Knudsen number flow) should then be 
scheduled for near-limb passes with a data accumulation rate high 
enough to detect possible density gradients suggested by Figs. 7, 

8, and 9. A suggestion by Barnes et al. (l97l) is possibly related 
to this discussion. These authors suggest that limb shocks can 
be caused by localized magnetic anomalies intrinsic to the moon. They 
further suggest that small length scales (geometric mean of electron 
and proton gyroradli) are ex]pected to promote instabilities which 
are, in turn, Induced by wave-particle interactions. Pertinent to 
the present study is the possibility that anomalous transport 
properties (i.e., "collisions") would be generated, 

EquaUy important would be detection of streamline patterns 
suggested by Pig. 10. This figure indicates that there could be a 
sudden deceleration along the subsolar streamllneB within a distance 
Of <-l radius from the surface. Scattering at higher energy densities 
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occuri3 near the surface in the direction of the on-coming flow, but 
the protons are quickly turned and deflected aroimd the limb» It is 
surprising to note that the dark side of the object is impacted by 
a very low, but finite, density of particles. The observations of 
Weugebauer et al. ( 1972 ) and Freeman ( 1972 .) noted above may, again, 
be relevant within the present context. 

We also point out that the present kinetic approach to the 
interaction of the solar wind with "email" obstacles (i.e., those 
celestial objects characteriaed by large Knudsen nuinbers) is suggested 
as an alternative to continuum and quasi-kinetic techniques. Indeed 
the results found in the Investlgatlcn suggest that, for a given 
proton thermal gyroradius, a shock wave would develop first at the 
limb, gradually "yiullding up - as we choose larger obstacles - to 
a *bonventional" detached bow shock as in Earth's case. Alternatively, 
for a given small obstacle , the effect would be the same if the 
solar wind were to get much hotter so that the gyroradius were to 
decrease gradually by several orders of magnitude. 

We state this important point in an alternate way, that is, by 
discussing an increase In the obstacle’s characteristic dimension 
In a slightly different way. The increase in size could certainly 
be gained by added material (l.e,, choosing a larger planet), but 
it could also be achieved by invoking the presence of an atmosphere, 
ionosphere, and/or a magnetosphere. The physical mechanism for 
new boundary conditions (ionopause, magnetopause, or even surface 
penetration) is a condition which could be Incorporated in related 
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studies. For example, surface absorption on‘the sunlit hemisphere, 
as in the case of the moon, could be incorporated in the expression for 
the number density of part ieles which are reflected from the 
surface (Eq. l4) . 

Finally, we should note that Wang (l97l) bas considered the 
magnetic field effect In his calculation of solar wind -moon inter- 
act ion; however, the condition of a completely absorbing surface 
is used. The reader is reminded that the magnetic field is not 
explicitly taken into account in the present study; although an 
estimate of its effect suggests that introduction of the Lorentz 
force would be confined to a skewing of the density contoux’s and 
an associated distortion of the flow pattern. Also, provision is 
not made in the present work for particle absorption (with no 
subsequent reemisslon)' in the sunlit hemisphere. 

Inasmuch as several statements have been made regarding the 
development of shock waves, "Mach cones^" and the like, provided 
the ctiaract eristic fluid scale length, is small relative to the 
obstacle’s scale length, a final comment is made to relate the 
present kinetic approach to its extreme "competitor"; continuum 
theory. It is well known that the "gain and loss" term in the 
kinetic description is given by the collision Integral. For free 
particle flow, this term is zero because there are no collisions 
of any kind. Similarly, in continuum theory - at least at the level 
of Euler’s approximation regarding the absence of transport prepei^ie 
such as viscosity - tlie collision integral is again found to be zerz 
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because the distribution function is found to be Maxwellian. There 
is, then, no contradiction in the - statement that kinetic theory, 
as discussed herein, preserves the characteristics of the continuum 
theory. 
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APPICWDIX 


We have employea the finite difference technique for obtaining 
nujjjerlcal rosults. We have chosen that "i" be the label for the 
r - coordinate and be the label for the 0 -coordinate, starting 
clockwise with Intervals of 15° in a range 0 f 6 ^ l 80 °. Also we 
have set up the origin of the coordinate system at the center of 
the body, therefore the local velocity distribution is expressed 

assumed that the local 

velocity distribution is Maxwellian forR 10 (which is effectively 
R = w) and 0 s 0 ^ iBO^, then ve compute each point successively 
by using Eq, (l) and Eq. (2). Newton’s equations of motion (i.e.. 
Equation (4) Is used here to track the particle trajectory 
due to the effects of the non-central force potential field. Because 
of the non-straight line characteristics, we can only fix our 
initial points on the boundary; the subsequent points, as said, 
were determined by Newton’s equations of motion. However, our 
grid points are fixed points, therefore, a Range-Kutta subroutine 
Is i^ised for the interpolation. 
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9 . 85 , Results are shown for an angularly-distributed surface 
electric potential. The surface repulsive electric field at 

the eubsolar point is taken to be 1 Vm~^. The result for no 
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7 . Constant density contours with non-central force field potential 
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9* Constant dencilty contours with central -force -fie Id potential 

for n - 5 euid S == 9*^5* 

10. Velocity vectors (or "streamline" pattern) with non-central- 
fpree fi._e Id potential for S - 9*^^5 in the vicinity of a 
"small" spherical ohstacle. Surface electric field at the 
subsolar point is 1 Vnn“-t, Note turning of the flow on the 
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CHAPTER V 

CONCLUSIONS AND RECOMMENDATIONS 


During fhis sfudy we invesHgated dynamicaS problems ranging from the 
earth's upper otmosphere to interplanetary space . The theories used to 
construct the models to seek the solutions for these problems also covered a wide 
range; from continuum theory of magnetohydrodynamics to discrete motion of 
particles, i kinetic theory of plasmos. 

It is not possible to obtain a complete solution due to the complexity of 
the mathemotical model of nature of the physical problems which we have 
attacked. Therefore, a number of appropriate assumptions are introduced. 

However, the results which we have obtained are reasonably accurate 
when compared with th© observations. There is still room for improvement of 
our models. Outlined oriefly ore possible future efforts. 

(1) The dynamical model for the responses of the neutral atmosphere 
is a very important subject in understanding the characteristics of the 
atmosphere from 90 Km to 200 Km. In the present study, we treat the 
atmosphere as a mixture of gases* The effects on each component of the 
constituent gases in the mixture is simply ignored. In order to understand 
the energy budget and concentration of the atmosphere, it is necessary 
to examine the dynamical equations for each constitment gas. This 
problem should be examined carefully. 

(2) Concerning the propagation of hydromagnetic waves in the ionosphere, 
only one particular case has been examined as discussed in Chapter III. 

Other important waves such as ULF, ELF and VLF need to be investigated. 

(3) Finally, the subject we have touched is the solar wind interoction 
with plonetary objects. We have neglected the interactions between 
ions and electrons, and ignored the wake and sheath structure. Therefore, 
if possible, these problems deserve particular attention in our future 
Spacelab program. 
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APPENDIX A 


List of Lectures Held on the Campus of The University of 
Alabama in Huntsville Resulting from this Contract 



P.O. Box 1247 
Huntsville, Alabama 35807 


Announcement of Lectuie 


THE RESTLESS GEOMAGNETIC FIELD 


Speaker: Dr. Wallace H. Campbell 

Time; I0:30a.m., Thursday, 20 December I973 

Place; The University of Alabama in Huntsville 

Research Institute Building, Room M-50 
Sparkman Drive, Huntsville, Alabama 


Abstract 


A review of the geomagnetic field changes occurring with periods between 0.5 sec 
and 1 month which reflect the dynamic behavior of our space environment. 


Biographical Information 

Currently Geophysicist, Geological Survey, Deportment of the Interior, Dr. Campbell 
received the B.S. (Louisiana State University, 1950), the M.A. {Vonderbilt University, 
1953), and the Ph.D. degrees (University of California at Los Angeles, 1959). He has 
taught at various institutions, lectured widely and published many scientific papers. 
Among his committee assignments have been the American Geophysical. Union, the 
International Association of Geomagnetism and Aeronomy, International Scientific 
Radio Union / and Inter-Union Commission on Solar Terrestrial Physics. Dr. Campbell 
has organized many of the national and international meetings he has attended as 
a delegate. 


For further informal Ion , call Dr. S. '! . V'/u at 895-6413. 



P.O. Box 1247 
Huntsville, Alabama 35807 


Announcement of Lecture 

THE REMOTE SENSING OF STRATOSPHERIC AEROSOLS FROM SATELLITES 


Speaker; 

Dr, Fred L. Bartman 

Time:' 

10:30a.m., Mondoy, 19 November 1973 

place: 

The University of Alobama in Huntsville 
Room M-50, Research Institute 

• 

Sparkman Drive, Huntsville, Alabama 


Abstract 


Members of the High Altitude Engineering Laboratory of the University of Michigan 
are studying the feasibility of determining stratospheric concentrations of minor 
constituents and pollutants from satellite measurements of infrared absorption during 
solar occultation. As part of this study the extinction of simplified models of 
stratospheric aerosols has been examined to determine: 1) How they might inter- 
fere with the measurements of molecular constituents, and 2) VVhether continuous 
monitoring of stratospheric aerosols may be corried out by satellite solar occulta- 
tion measurements. The results of this study of aerosols ore discussed after the 
aerosol model and its optical characteristics are described in detail. 


Biographical Sketch 


Dr. Bortman received the Ph.D. from the University of Michigan, where be is 
currently Director of the High Altitude Engineering Laboratory and Professor in 
the Departments of Aerospace Engineering and Atmospheric and Oceanic Science. 
An expert in high altitude atmospheric sciences, he has published numerous 
papers and served as o consultant to many different governmental agencies. 


For further informal Ion, call Dr. 5. T. Wu, 395-6413. 



P.o. Box 1247 
Huntsville, Alabama 35007 

Announcemo nt of Lecture 

IONOSPHERIC STORMS IN LOW LATITUDES AND 

INTERPLANETARY SECTOR FIELD EFFECTS ON THE POLAR CAP 

Speaker: Dr. S. Mahushlto 

* 10:30 a. m., Friday, 9 February 1973 

The University of Alabama in Huntsville 
Research Institute, Room M-50 
Sparkman Drive 
Huntsville AL 

Abstract 

The ionization density increase in the F2 layer, on the overage, over the magnetic 
equator during geomognetic storms has long been a puzzle. This Is explained by a 
decrease in the fountoin effect mainly caused by a change of the E-W electrostatic 
field (hence the vertical plasma drift). 

It has recently been reported by several scientists that the interplanetary magnetic 
sector fields causa the polar cap geomagnetic field changes. This phenomenon is 
re-examined In detail and is explained by a shift of the Sq current pattern. 

Biographi cal Information 

Currently Senior Scientist, High Altitude Observatory, National Center for Atmospheric 
Research, Professor of Astro -Geophysics, Univ. of Colorado, and consultant to NCAA, 
Dr. Matsushita holds the M.Sc. and the Dp. Sc. (Kyoto Univ.). He was a lecturer at 
Kyoto Univ. and o member of the research staff, University College, London, before 
joining the staff of the Univ. of Colorado in 1955. He is the recipient of the Tanakadate 
Award, Society of Terrestrial Magnetism and Electricity, 1950, and the Scientist Award, 
Science Research Society of America, 1963. Dr. Matsushita is co-author of two books, 
has boon editor of journals and proceedings, and published over 120 papers. Principol 
research has been on solar terrestrial relationship and physicol behavior ionosphere, 
magnetosphere and the earth’s magnetic field. He discovered lunar effects on lower 
ionosphere and the* special ionospheric zone of sporadic E over the magnetic equator. 


For further information, please coll Dr. S. T. Wu - 895-6413 
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Announcement of Lecture 


SOME RECENT ADVANCES IN TIME SERIES ANALYSIS 


Speaker; 

Dr. Emanuel Pa rzen 

' 

Chairman, Department of Statistics 
State University of New York at Buffalo 

Time: 

10:30a.m., Friday, 3 November 1972 

Place: 

The University of Alabama in Huntsville 
Research Institute, Room M-50 
301 Sparkman Drive, N.W. 

Huntsville, Alaboma 


Abstract 

This survey of current research will include a discussion of: estimation of 
covariances and inverse covariances of stationary normal time series; 
autoregressive spectral estimotlon; efficient closed form estimation of the 
parameters of moving average schemes and mixed autoregressive moving 
average schemes; and estimation of parameters of models of signal plus noise. 

Biographical In formation 

Dr. Parzen received the B.A. (1949) from Horvard, the M.A. (1951) and the 
Ph.D. (1953) degrees from The University of California at Berkeley. Among the 
positions he has held are Research Scientist, Hudson Labs., Columbia University; 
Guest Professor of Mathematics, Imperial College, London; Visiting Professor of 
Statistics, MIT; Professor of Statistics, Stanford University; and a Fellow, ISM 
Systems Research Institute, New York City, before becoming Professor and 
Chairman, Department of Statistics, State University of New York at Buffalo. 

Dr. Parz.en holds membership In numerous associations and has been referee for 
many journals. Associate Editor, SIAM Control Journal and Editor, Time Series 
Analysis Series of Holden Day, he has authored over 35 papers and, among others, 
the following books: Modern Probabilit y Theory and Its AppI ications, Stochastic 
Processes and Time Series Analysis Papers. 


For further Information, call: Dr. S. T. Wu, 895-6413 
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Announcement of Lecture 

DISSIPATION OF HYDROMAGNETIC WAVES 


Speaker: 


Dr. R. J. Hung 


jlnie: 10:30 a.m., Monday, 19 June 1972 

Place: The University of Alabama In Huntsville 

Research Institute Conference Room 
Sparkman Drive 
Huntsville, Alabama 35807 


Abstract 


Propagation of hydromagnetic waves is modified by transport phenomena due to 
Coulomb collisions. CoKisIonal effects can vary from region to region in a plasma 
of interest. In particular, woves of a given period may see one region of a 
plasma as collisionless (in the sense that the wave period Is short compared with the 
Coulomb collision time) and another region of the plasma as collision-dominated. 
For exomple, a tronsition from collisional to collisionless behavior may take place 
for a wave propagating upward in the ionosphere or solar corona, because the 
Coulomb collision time decreases with altitude. It is therfore of considerable in- 
terest to investigate how hydromognetic waves behave In the whole plasma. In- 
cluding the collisional, intermediate and coliislonless regions, and in particular, 
how the waves are dissipated. The work reported here is concerned with the dis- 
sipation of hydromagnetic waves, especially In the transitional outer solar corona. 
The theory Is equally applicable to the Ionosphere and other areas of space 
physical and astrophysica! problems. 


Biographical Information 

Dr. Hung received the B.S. degree In 1957 (Notional Taiwan University), the 
M.S. in 1966 (University of Osaka) the Ph.D. in 1970 (The University of 
Michigan). A National Research Council Postdoctoral research associate, he 
is currently working in the TheoreMcai Studies Branch of Ames Research Center. 
With his background in plasma physics mognetehydrodynamlcs and fluid 
mechanics, he has authored various papers and received several honors and 
awards. 

For further ii^formatlon, please call Dr. S.T. Wu - 895-6413 
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RESEARCH INSTITUTE, 
DIVISION 01- GHADUAte 
PROGRAMS AND, RESEARCH 

P.O. Box 1247 
Huntsville, Alabama 35807 


Announcement of Lecture 
STUDIES OF ATOMIC COLLISIONS IN THE LABORATORY 


AND IN THE IONOSPHERE 
Speaker: Dr. M. A. Biondi 

Xime: 10:30 a.m., Friday, 19 May 1972 


The University of Alabama in Huntsville 
Research Institute, Room M-50 
Sparkman Drive 
Huntsville, Alabama 35807 


Abstract 

The path connecting laboratory studies of relevant atomic collision processes 
with models of ionospheric behavior is illustrated using dissociative electron- 
ion recombination and charge transfer as examples. Microwave and optical 
interferometer techniques are applied In laboratory studies to determine the 
nature of the electron-ion recombination process and its dependence on 
electron energy. Drift tube - mass spectrometer techniques are used in the 
study of low energy ion-molecule reactions such as non-resonant charge transfer. 
The results of these laboratory studies ore applied to a model of ionospheric 
behavior under normal and artificially perturbed conditions. The predictions of 
the model are compared with Interferometric studies of 6300A 01 nightglow 
emission line profiles ond with photometric studies of 6300A Intensity changes 
produced by powerful ground-based transmitters. 


Biographical Information 

Dr. Biondi received the Ph. D. In physics in 1949 from the Massachusetts 
Institute of Technology. From 1949-52 he was resident physicist, later, Manager 
of the Physics Department, Westlnghouso Research Laboratory. Currently 
Professor of Physics at the University of Pittsburgh, he has been consultant to 
various agencies. Including NSF and NRL, and published numerous papers. 

Areas of Interest are atomic physics, interactions and reactions, involving electrons, 
ions and excited atoms, solid state physics, electromagnetic properties of metals 
at liquid helium temperatures. 


For further information, please call Dr. S. T, Wu - 895-6413 
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Announcement Lecture 


INTERACTION OF ELECTRONS WITH 
WATER VAPOR WITH ATMOSPHERIC APPLICATIONS 


Speaker: Dr» John J, Olivero, Jr* 

Postdoctoral Research Associate-Instructor 
(Physics and Astronomy) 

University of Florida 

Time; 10;30a.m., Friday, April 14, 1972 

Place: The University of Alabama in Huntsville 

Research Institute Conference Room 
Sparkman Drive 
Huntsville, Alabama 


Abstract 


We examine the consequences of electron impact on water vapor in terms of the 
microscopic details of excitation, dissociation, ionization, and combinations of 
these processes. Basic electron impact cross section data are assembled in many 
forms and incorporated into semi analytic functions suitable for onolysis with digital 
computers. Energy deposition in water vapor is discussed and the energy loss func- 
tion (“i dE/dx) is presented along with the "eV per ion pair" and the efficiencies 
of energy loss in various processes. We consider several applications of electron- 
water vapor interactions in the atmospheric sciences. We discuss H 2 O comets, 
aurora and airglow, and lightning. In general we find that other processes are the 
more likely source of H*, O*, and OH* emissions with the probable exception 
of lightning. 


Biogrgphicol Information 


Currently Postdoctoral Research Associate-Instructor (Physics and Astronomy) , 
University of Florida, Dr. Olivero holds the M. S. degree from College of William 
and Mary and received his Ph.D. from University of Michigan in 19701 He has 
been employed by NASA-Langley Research Center in Virginia^nd the University 
of Michigan. Dr, Olivero is the author of several publications. 


For further Information, call Dr. S, T. Wu ot 895-6413, 



RESEARCH INSTITUTE. 
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PROGRAMS AND RESEARCH 

P.a Box 1247 
Huntsville, Alabama 35807 


Announcement" of Lecture 


THE NEGATIVE BINOMIAL DISTRIBUTION: DIFFERENT POINTS OF VIEW 


Speaker: 

Dr. G. P. Patil 

Time; 

10:30 a.m., Monday, 27 March 1972 

Ploce: 

The University of Alabama in Huntsville 
Research Institute, Room M-50 
Sparkman Drive 
Huntsville, Alabama 35807 


Abstract 


The negotive binomial distribution has been studied extensively in the post, and it 
has been found that data arising in various situations in many fields fits this 
distribution ♦ Naturally, then, one wonders about the underlying chance mechanisms, 
and it turns out that there are several such mechanisms which have no apparent 
relationship. The emphasis of this paper is on bringing together in a single notation 
stochastic models generating the negative binomial distribution (NCBD), and the 
basic cboracter of the paper is that of a systematic review and exposition of the 
Qvoilable results on the subject. A few new results and comments given will, 
hopefully, clarify the structure of this important distribution. With the exception 
of the important o~truncated negative binomial distribution, distributions formed 
by randomly stopped sums, mixtures and other modifications of the negotive 
binomial distribution are not discussed.' 

Biogrophicql Information ' 

Dr. Patil received the M.Sc. (1955) from Poona Lhiversity and the M.S, and Ph.D. 
(1959) from the University of Michigan. He has served on the faculties of the 
universities of Michigan, McGill and Wisconsin and presently is Professor of 
Mathematical Statistics at Pennsylvania State Lhiversity. He has been consultant 
with USDA, Forest Service, NIH, HEW, and the Laboratory of Applied Studies. 

He is the author or co-author of "Dictionary and Bibliography of Discrete Distributions," 
some 50 research publicotlons ond editor of "Random Counts in Scientific Work" 
and"Stol istical Eculogy" and director of severed intcinotional symposia. 
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For further information, please call Dr. S. T. Wu - 895-6413 
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Huntsville, Alabar^a 35807 

Announcement of Lecture 

MEASURING WAVES, ESTIMATING WAVE SPECTRA AND SPECIFYING AND FORECASTING WAVES 

GIVEN THE WIND FIELD AND WIND FIELD FORECASTS 


Speaker: 

Dr. Willord J. Pierson 

Time: 

' 10:30 a. m., Thursday, 9 March 1972 

Place; 

The University of Alabama in Huntsvill 


Research Institute, Room M-50 


Sparkman Drive 


Huntsville, Alabama 


Abstract 


Waves have been measured by the Tucker Shipborne wave recorder ond by stereo- 
photogramme trie, buoy and airborne laser techniques. From these measurements, 
the spectra of the gravity wave components for so-called "fully developed seas," 
for growing seas as a function of fetch and duration, and for dead seas and swell 
have been found. From these data, and from theories on wave generation, non- 
linear interactions, wave dissipation and wave propagation, a computer based 
procedure for specifying and forecasting waves on the open ocean in terms of 
their spectra, given wind fields ond wind field forecasts, has been developed. 
Such computed spectra can be used to determine the wave climate for a given 
area of the ocean and to predict how ships, towed objects, such as spent rockets 
cyclinders, and buoys would behove in the waves. 


TRE UNIVERSITY OF ALABAMA IN HUNTSVILLE 



Biographical In formation 

Professor Pierson received the B.S. degree from the University of Chicago in 1944 
and the Ph.D. degree from New York University in 1949, when he carno o mem- 
ber of the faculty; he is Professor of Oceanography currently. Editor of the 
Journal of Atmospheric Science, he is a member of the American Geophysics 
Uhion (vice president, section on oceanography 1965-57), American Meteorology 
Society, American Astronoutlcal Society, Marine Technology Society and Sigma 
Xi. He has published two books and numerous articles. Principal research jn- 
torosh include motion of ships in waves, v/avo analysis by time series, ond fore- 
casting by computer tocliniques, use of radar on spococraft to study oceans, 
turbulence, ocean currents. 


For further informarlon, please coll Dr. S, T, Wu - 895-6413 
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PROGRAMS AND RESEARCH 
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Hunc&vtile, Alabama 35807 


Announcement of Lecture 


DAILY VARIATION OF THE IONOSPHERIC F2 EQUATORIAL 
ANOMALY IN THE AMERICAN AND ASIAN SECTORS 


Speaker: Dr, David N. Anderson 

CIRES Visiting Fellow 
University of Colorado 

Time: 10:30 a, m., Thursday, December 16, 1971 

Place: The University of Alabama in Huntsville 

Research Institute, Room M-50 
Sparkman Drive 
Huntsville, Alabama 


Abstract 

When observed noontime values of the maximum electron density in the ionospheric 
F2 region are plotted as a function of magnetic latitude, a curve is produced which 
has two peaks, one on either side of the dip equator at + 16° dip latitude. The 
present work mainly will discuss the theoretical investigations of the daily variation 
of this latitude distribution in NMAX(F2) (the so-called Appleton cr equatorial 
anomaly) and specifically attempts to account for the longitudinal differences ob- 
served between the American and Asian sectors. Other investigations include the 
effects of the neutral wind in produciiYg north-south asymmetry in the anomaly ond 
the solar cycle effect on the anomaly's daily variation. 

Biographical Information 

Currently CIRES Visiting Fellow at the University of Colorado, Boulder, Dr. Anderson 
holds the M.S. degree from North Carolina Stote University and received his Ph.D, 
from University of Colorado in 1971 . He has been employed by Los Alamos Scientific 
Laboratory In Now Mexico, NCAR and Dept, of Commerce in Boulder, Colorado. 
Member of American Geophysical Union, Phi Koppa Phi, NDEA Fellowship ('62-'64) 
North Carolina State, HAO Graduate Assistantship, University of Colorado. He 
has published some technical jxjpers. 


For further information, call Dr, S, T. Wu at 895-6413 
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P.O. Box 1247 
Huntsville. Alabama 35807 


Announcement of Lecture 


MAGNETISM OF THE EARTH IN SPACE 

Speaker: Dr. S. Matsushita 

Time; 2 p.m., Thursday, 18 November 1971 

piece: The University of Alabama in Huntsville 

Research Institute, Room M-50 
Sparkman Drive 
Huntsville, Alabama 

Abstract 


THE UNIVERSITY OF ALABAMA IN HUNTSVILLE 



Behavior of the earth's Ionosphere and magnetosphere estimated from the viewpoint 
of geomagnetic variation fields are reviewed. Key observations needed to resolve 
the remaining theoretical questions concerning magnetically quiet conditions are 
discussed. Concerning magnetic disturbances, the current state-of-the-art and a 
new idea based on recently available data are presented, including suggestions to 
find missing links. 


Biographical Information 


Currently Senior Scientist, High Altitude Observatory, Notional Center for Atmospheric 
^.esearch. Professor of Astro-Geophysics, Univ. of Colorado, and consultant to NOAA, 
Dr. Matsushita holds the M.Sc and the Dr. Sc. (Kyoto Univ.). He was a lecturer at Kyoto 
Univ. and a member of the research staff, University College, London^before joining the 
staff of the Univ. of Colorado in 1955. He is the recipient of the Tanakadate Award, 
Society of Terrestrial Magnetism and Electricity, 1950, and Hie Scientist Award, Science 
research Society of America, 1963. Dr. Matsushita is co-author of two books, has 
ceen editor of journals and proceedings, and published over 120 papers. Principal 
.*eseorch has been on solar terrestrial relationship and physical behavior ionosphere, 
ir-jagnetosphere and the earth's magnetic field. He discovered lunar effects on lower 
’onosphere and the special ionospheric zone of sporadic E over the magnetic equator. 


'O' further information, please call Dr, S. T. Wu - 895-6413 
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Announcement of Lecture 


DYNAMIC EFFECT OF SUDDEN HEATING IN THE AURORAL ATOMOSPHERE 


Speaker 

Dr, Paul B. Hays 

Time: 

lOa.m,, Friday, 22 October (97I 

Place; 

The University of Alabama in Huntsville 
Graduote Studies Building, Room llO-lll 
Sparkman Drive 
Huntsville, Alabama 


Ab$troct 

The auroral atmosphere is, during times of magnetic storms, subjected to intense 
local heating. This energy addition results both from particle impact and electric 
fields; however, the electric field heating appears to be dominant. The neutral 
atmosphere responds to local heating by setting up cellular dynamic systems which 
cause changes in thermospheric composition. These composition changes have 
been detected and are compared with theoretical predictions. 


Biographical Information 

Currently Professor of Aerospace Engineering at the University of Michigan, from 
which he received the Ph.D, degree, Dr. Hays did post doctoral work at Queens 
University, Belfast, A leading authority iri atmospheric science. Dr, Hays has 
published numerous papers and books In his Field. 


For further information: Call Dr. S. T. Wu - 895^6000 
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Annooncoment of Lecfuro 


SPECTRAL MOMENT ESTIMATION 


Spoalcer: 

Tlmo: 

place: 


Dr. Kennerh S. Miller 

I0a*m., Thursday, 13 May 1971 

The University of Alabama in Huntsville 
Research Institute Lecture Room, M-50 
301 Sparkman Drive 
Huntsville, Alabama 

Abstract 


The problem to be discussed Is estimating the moments of the spectral density of a complex 
stationary process. This problem is of interest in such diverse fields as tropospheric com*- 
munication channels, turbulent velocity phenomena (ionized wakes) and the estimation of 
arterial blood velocities using Doppler techniques. 

Some of the details have appeared in published papers. This brief lecture will outline the 
general philosophy governing the choice of suitoble estimators, explore some of the ramifi- 
cations of the problem, and indicate how to compute the statistics of various estimators. 

Biographical Information 

Dr, Kenneth S, Miller, a native Now Yorkel*, received the B,S, degree in chemical 
engineering and the M,A. and Ph.D. degrees in mathematics from Columbia University. 
Post-doctoral work was don© at th® Institute for Advanced Study, Princeton, in 1950. 

A radar officer in the U. S. Navy during World War U, Dr. Miller was, successively, 
instructor, assistant professor, associate professor and Professor of Mathematics at New York 
University, During this time h® also acted as a consultant to various industriol and 
governmental agencies on problems associated with system analysts, noise, mathematical 
machines, ond applied mathematics. He is the author or co— author of over fifty research 
papers and fifteen books. Currently he is at the Riverside Research Institute, New York, 
New York, and as Adjunct Professor of Mathematics at Fordhom University. 

Dr. Miller is a member of Sigmo Xi, Tau Beta PI, Pi Mu Epsilon, the American Mathematical 
Society, and a senior member of th© IEEE. 


For further information: Call Dr. S. T. Wu - 895-6413 
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Announcement of Lecture 

I III I III '■ ■ I ■ I ■ I ^ .1 


CHEMISTRY OF THE UPPER ATMOSPHERE 
AND THE ROLE OF TRANSPORT PROCESSES 

Dr. Richard A. Craig 

10 a. m., Tuesday, A/brch 9, 1971 

The University of Alabama in Huntsville 
Research Institute Lecture Room, M-50 
301 Sparkman D rive 
Huntsville, Alabama 


Speaker: 

Time; 

Place: 


Abstract 

The important chemical reaction affecting the neutral composition of the upper stratosphere, 
mesosphere, and lower thermosphere will be discussed. The role of atmospheric transport 
processes will be emphasized and the problem of taking these into account elucidated. 

Biographical Information 

Dr. Craig received the A.B. degree In astronomy from Harvard in 1942, the M.S. ond Sc . D. 
both In meteorology from M.I.T., in 1944 and 1948 respectively. Currently Professor and 
Choirmon, Department of Meteorology at Florida State University, Dr. Craig served as a 
research assaxiatte at Woods Hole Oceanographic Institution, a research fellow at Harvard 
College Observotopy, and a researeh meteorologist and bronch chief at the Air Force Cam- 
bridge Research Center. Active in professionol societies. Dr. Craig has been associate 
editor or the Journal of Meteorology ond the Journal of Geophysical Research. He has 
been a consultant to the Notional Science Foundation, White Sands Missile Ronge ond NASA. 
He has published widely In various journals, co-authored a number of books and outhored 
The Edge of Space; Exploring the Upper Atmosphere (Doubledoy). 


For further information: Call Dr. S. T. Wu 895-6413 
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THE RESEARCH INSTITUTE 
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Huntsville, Alabama 35007 


Announcement of Lecture 


SOME ASPECTS OF ATMOSPHERIC WAVES IN REALISTIC ATMOSPHERES 


Speaker: 

Time: 

Place: 


Dr. Richard Seigmund Lindzen 

10 a, m./ Tuesday, 2 February 1971 

The University of Alabama in Huntsville 
Research Institute Lecture Room M-50 
301 Sparkman Drive, N.W. 

Huntsville, Alabama 


Abstrac t 

The theoretical behavior of tides and gravity woves will be discussed, accounting for the 
effects of viscosity, thermal conductivity ond mean shears. Special emphasis will be given 
to the semi-individual tide and to atmospheric free oscillations. 

Biographical Information 


Dr. Lindzen received the B.A, degree in physics (I960), the M.S. (1961) and the Ph.D. (1964) 
in applied mathematics from Harvard. The following year he wos appointed a research associate 
at the University of Washington and, in 1965, a NATO post doctoral fellow at the Institute for 
Theoretical Meteorology at the University of Oslo. He received appointments at the National 
Center for Atmospheric Research as a research scientist from 1966 to the end of 1967, returning as 
a lecturer the following summer. He hos been a visiting lecturer at the University of California 
at Los Angeles, a visiting professor at Tel Aviv University, and Associate Professor of Meteorology 
ot the University of Chicago since 1968. 

A member of vorious organizotlons ond committees, Dr. Lindren is the recipient of the N.C.A.R. 
Outstanding Publicotion Award, the A. M.S. Meisingcr and A.G.U. Macelwane Awards and an 
Alfred P, Sloan Fellowship. Ho has published extensively, both in the United States and in 
foreign journals. 


per further informotion- Coll Dr. S. T. Wu, 895-6413 
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Announcement of Lecture 


THE STABLE AURORAL RED ARC AND ITS 
INTERACTION WITH THE NEUTRAL ATMOSPHERE 

« 

Speaker: Dr. R. G. Roble 

Time: 10 a, m., Tuesday, December I, 1970 

P^ace: The University of Alabama in Huntsville 

Research Institute Lecture Room M-50 
301 Sparkman Drive, N.W. 

Huntsville, Alabama 

Abstract 

The stable auroral red arc which occurs in subauroral latitudes during geomagnetic storms is 
examined in the light of recent observations. The various excitation mechanisms which have 
been postulated to explain the red arc are examined critically. Satellite observations of 
electron temperature and the topside electron density structure are used to calculate a red 
arc considering thermal conduction in the electron gas as the only energy source. The calcu- 
lated A. 6300 emission rate, geographical position, and horizontal exten are compared with 
photometric data obtained from airglow observatories and the results show good agreement. 

The calculated neutral gas temperature response to electron heating within the arc is shown 
to be small in agreement with A 6300 doppler temperature measurements made with a Fabry- 
Perot interferometer on the red arc of October 3)/November I, 1968. 

Biographical Information 

Raymond G. Roble received B.S.E. degrees in Engineering Physics and Engineering Mathematics 
(1957), the M.S.E. degree in Mechonical Engineering (1961), and the Ph.D. In Aeronomy (1969), 
University of Michigan. After serving three years with the U. $. Navy and later working with 
Bendix Research Laboratories, Dr. Roble returned to the University of Michigan High Altitude 
Engineering Loboratory. He designed and built the Airglow Observatory, which has a 6** diameter 
Fabry-Perot Interferometer ond an airglow photometer, used fo study the stable auroral red arc. 

• Currently Dr. Roble Is an atmosphe:ic scientist ot the National Center for Atmospheric Research, 
Boulder, Colorado. 


For further information: Coll Dr. S, T. Wu, 895-6000 
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Huntsville. Alabama 35807 


Announcemen^ of Lecture 


"THE COMPOSITION OF UPPER ATMOSPHERE" 


Speoker. Professor Alfred Nier 

10 a. m., Tuesday, November 13, 1970* 

Place: tl 1 1 • . 

Ihe University of Alobama In Huntsville 

Research Institute Lecture Room M-50 
301 Sparkman Drive, N.W. 

Huntsville, Alabama 

Abstract 

spectrometers, carried on sounding rockets and satellites, hove been soccessfollv emo/oved 
or studying the neutral and ionic composition of the earth's atmosphere obout 100 km The ^ 
meosurcment problem is complicated by the high speed of the vehicles, but methods have been 
eloped for minimizing the difficulties. Above 100 km neutral N2, 02,0 and Ar ore in dif- 

loi: »• — 

Biographical Info rmation 

InTic" r'Lungln'^i^^^^^^^^^ Th' D B-A- (1931) and M. S. (1933) 

rently Dr. NIer I, Regents' Profc»o‘r of Phy, ics ofX' uTver i.y ex- 

cessively wos Assistant Professor (1938-40), Associate Professor a940-47)"Cf;sror1ri^"o794'lt 

r 

,,, Ke.ie;corporotiorrpoV:;7hr;M:;h^^^^ 


'Originally announced date: 23 October 1970 


For further information; Coll Dr. S. T. Wu, 89.'j-64l:3 
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Annouhcement of Lecture 


"STRUCTURE AND VARIATIONS OF THE HETEROSPHERE" 


Speaker: Dr. Luigi G. Jacchia 

Time: 10 a . m. , 29 October 1970 


Place: 


The University of Alabama in Huntsville 
Research Institute Lecture Room M-50 
301 Sparkman Drive N.W. 

Huntsville, Alabama 


Abstract 

Atmospheric composition is kept nearly homogeneous by mixing to o height of about 90 km (the 
homosphere); above 100 km we have diffusive separation of the atmospheric constituents (the 
heferosphere) . The heterosphere is dominated by two energy sources: solar EUV and the solar v/ind; 
the manyfold, large variations that are observed in this region can all be related, directly or in- 
directly, to variations in the Intrinsic intensity of these two sources and to variations in the loca- 
tion of the sources relative to the observer. Atmospheric models can be successfully constructed, 
in which the temperature and density variarions are empirically reloted to solar and geomagnetic 
parameters. Some of the more recent findings concerning the structure and variations in the hetero- 
sphere will be presented. 

Biographical Information 

Born in Trieste, Italy, 1910, Dr. Jacchio received the Ph.D. degree from the University of Bologna 
in (932. Currently a Physicist at the Smithsonian Astrophyslcal Observatory, he wos a D.I.C. stafi 
member (1943-49) and a Research Associate in Electrical Engineering (1949-53), at M.I.T. From 
1941-45, Dr. Jacchia did defense work, foreign language broadcasting and monitoring and was a 
Scientific Consultont to the Office of War Information. Previously, he was a Reseorch Associate 
at Harvard College Observatory, held o fellowship of the Society for the Protection of Science 
and Learning for Research, University of London Observatory, was a lecturer. University of Bologna, 
and an astronomer, University of Bologna Observotory, 


For further information: Call Dr. S. T. Wu, 895-6413. 



THE UNIVERSITY OF ALABAMA IN HUNTSVILLE 



TUt HlSLAFICmNSIITUTf 

P.O Box 1947 
Huntevtiio. Aiabama 35B07 


Atinounccnienf of Lecfure 


••UPPER ATMOSPHERIC SPECTROSCOPY" 


Speaker: 


Dr. A. E. S. Green 


Time: 10 a. m,, Thursday, October 15, 1970 

Place: The University of Alabama in Huntsville 

Research Institute Lecture Room M-50 

, 301 Sparkman Drive, N.W. 

Huhtsville, Alabama 


Abstract 

Spectroscopy and atomic physics are so intimately entwined with planetary aeronomy 
that progress in this geophysical-astronomical discipline in large measure is dependent 
upon progress in certain phases of those basic disciplines. In this lecture, we illustrate 
the interdependence by a discussion of photon, electron and proton stimulation of dayglow 
and auroral spectra. 


Biographical Informotion 


Dr, Green received the B.S. from the City College of New York (1940), the M.S. 
from the Californio Institute of Technology (1941), and the Ph.D. degree from the Univer- 
sity of Cincinnati (1948), all in physics. Manager of Space Science Laboratory, General 
Dynamics (1959-63), Dr. Green has held academic apfjointments at the University of Cin- 
cinnati as Associate Professor (until 1953), at Florida State University as Professor and 
Scientific Director, Tandem Van deGraoff Progrom (until 1959), and at the University of 
Florida as Graduate Research Professor (since 1963); ho served os Operations Analyst and 
Physicist with the U. S. Air Force (1944-45). Dr. Green held various visiting appoint- 
ments including Oak Ridge Notional Laboratory, Los Alamos, Aerospoce Corporation, 
and Stanford Linear Accelerator. 


For further information: Coll Dr. S. T. Wu, 895—6413. 
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“UPPER ATMOSPHERIC RESPONSE TO TRANSIENT HEATING" 


Speaker; 


Dr. Gary E. Thomas 


Time; 


10;30A.M., Friday, September 25, 1970 


place: Research Institute, The University of Alabama 

" In Huntsville, Lecture Room M-50, 301 

Sparkman Drive, N. W., Huntsville, Ala. 


Abstract 

The time -dependent response of the upper atmosphere to transient heat sources is 
considered. The basic problem is that of heating of a compressible, heat-conducting 
fluid, which is described in the one -dimensional case by an analytic solution. Com- 
parisons with satellite drag data of such first-order solutions are shown to be useful in 
determining energy requirements and in detennining some constraints on the spatial dis- 
tribution of the heating. Recent OGO-6 mass spectrometer and interferometric tempera- 
ture measurements show that otmospheric disturbances during geomagnetic storms are 
much more prominent in the ourorol zones. These results suggest that joule dissipation 
of auroral currents are important contributors to orbital perturbations of satellites. 


Biogrophicol Infor mation 

Dr. Thomas was born October 25, 1934, in Lookout, West Virginia. He received 
the B. S. in Physics at New Mexico State University in 1957 ond the Pb.D, in Physics 
at the University of Pittsburgh in 1963. He received a one year post-doctoral appoint- 
ment at Service d'Aeronomie, C.N.R.S., Paris, France 1962-63, He served as IstLt., 
Sig. C., U. 5. Army Electronics Laboratory, Ft. Monmouth, New Jersey, 1963-64. 

He was a Research Scientist at Los Alamos Scientific Laboratory, Los Alamos, New Mexico, 
in 1964; a member of the technical staff, Space Physics Loboratory, Aerospace Corporation, 
El Segundo, Collfornio from 1965-67; Assistant Professor (1967-69) and Associote Professor 
(1969 -present), Astro -Geophysics Department, University of Colorado, Boulder, Colorado. 
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EFFECTS OF DISTURBANCE ELECTRIC CU*RRENT ON THE UPPER ATMOSPHERE 

Dr, Keifh D. Cole 

10:30 A,M,, Tuesday, May 26, 1970 

Research Instll'ul'e, Room M-50 
The University of Alabama In Huntsville 
Huntsville, Alabama 

Abstract 

Sotr.o likely sources of electric currents and fields In the auroral ionosphere oro 
reviewed. The dissipation of these currents and fields is then developed. Among otnor 
things, the following two major effects aro shown to occur: (i) Increases of scale height 
end Vomperaturos at the site of the auroral electrojet at altitudes above about 100 km. 
and (i i) Generation of high wind speeds i.n the auroral atmosphere. 



Dr, Keith D. Cole was bom in Australia and received the Doctor of Science degree 
from the University of Queensland. In 1956 ha was a member of the Antarctic Expedition 
to Maccuorie Island. From 1957-63 ho was Theoretical Physicist, Australia National 
Antartic Research Expeditions, He has held the position of Research Associote with both 
the Universii-y of Chicago and the University of Colorado, Dr. Cole is presently a NAS-NRC 
Senior Post-Doctorcl Follow at Goddard Space Flight Center, Greenbelt, Maryland, on 
leave from La Trobc University, Victoria, Australia, where he Is Foundation Professor of 
physics. He has published over 50 papers in international journals, most on theory of mag- 
netic dlsturbanco, 

* it it it 

There is no charge for admission to the lecture and no pro-registration is required , For 
cdditior.al information, please coll Dr, S. T. Wu, Associate Professor of Engineering, The 
University of Alabama in Huntsville, Telephone No. 895-6177, 


Speaker; 

Time: 

Place; 
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THE FLARE FORECASTING FIASCO - IS EVERYBODY LOOKING AT THE SAME SUN? 


Speaker: 


Dr. Frederick W, Ward, Jr. 


Time: 


10:30 A.M., Tuesday, May 19, 1970 


p\ace: Science & Engineering Building, Room 127 

The University of Alabama in Huntsville 
Huntsville, Alabama 


Abstract 


The passionate dreom of the oppressed solar flare forecoster is the scientific discovery 
of the cause, the energy source and the observable antecedent conditions for the solar flare. 
This "tunnel vision" has had a profound and deleterious effect on the state-of-the-art. It 
has diverted attention from mony promising approaches and relegated them to limbo. More 
importantly, it has actually inhibited the scientific search for the underlying physical mech- 
anisms. Examples of neglected aspects of the problem cover the entire range from observing 
techniques and equipment to the ovailability and accuracy of the archived scientific data. 

A cursory inspection of this data leads to an obvious conclusion - there are at least two suns 
up there. 

Biographical Informotion 


Dr. Frederick W. Ward, Jr., of the Space Physics Laboratory, Air Force Cambridge 
Research Laboratories, Bedford, Massachusetts, , received the B.S., M.S., and Ph.D. degrees 
in meteorology from M.l.T. He directed the AFCRL Space Forecasting Program from 1966-69. 
His main research is on the problem of the general circulation of the solar atmosphere and the 
maintenance of the differential rotation. 
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List of Publications Resulting from this Contract 


(1) "Propagation of Hydromagnetic Waves in the Upper F~2 Region/ " 
Planet. Space Sci., Vol . 21, 1973. 

(2) "Kinetic Theory Analysis of Solar Wind Interaction With Planetary 
Objects, " in Photon and Particle Interaction With Surfaces in Space 
(Ed. by R. J. L. Grard), D. Reidel Publishing Company, Dordrecht- 
Holland, 1973. 

(3) "The Dynamical Responses of the Thermosphere Due to a Geomagnetic 
Storm," AIAA Paper No. 74-217. 

(4) "An Analysis of the Upper Atmosphere Wind Observed by LOGACS, " 
planet. Space Sci., 1974. 

(5) "Joule Heating and Winds Due to Geomagnetic Disturboncs, " J. 
of Atmospheric and Terrestrial Phy., Vol. 36, 1974. 



